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ABSTRACT 


A  brief  review  of  the  classical  theories  of  rigid  perfectly 
plastic  and  elastic-perfectly  plastic  conti nua  is  presented  and  the 
invariance  requirements  of  constitutive  equations,  which  lead  to  the 
kinematically  correct  form  of  the  Prandtl -Reuss  equations  are  discussed. 

The  kinematics  and  thermodynamics  of  elastic-plastic  de¬ 
formation  of  materials  which  may  undergo  finite  elastic  strains  are 
discussed  and  a  plastic  flow  rule  obeying  the  necessary  invariance 
requirements  is  proposed  for  these  materials.  Relevant  topics  in  the 

theory  of  finite  strain  elasticity  are  reviewed. 

Two  intepretations  of  the  von  Mises  yield  condition  are 
given  for  a  material  which  may  undergo  finite  elastic  strains,  and 
the  associated  plastic  flow  rules  are  found. 

Numerical  solutions  are  obtained  for  the  stress  deviator 
tensor  in  a  cuboid  undergoing  simple  shear.  The  two  elastic-perfectly 
plastic  materials  considered  exhibit  neo-Hookean  elastic  behavior  prior 
to  yielding  and  during  unloading.  Both  interpretations  of  the  von 
Mises  yield  condition  and  the  proposed  plastic  flow  rule  are  used. 

The  numerical  results  found  for  the  simple  shear  problem 

O' 

are  used  to  obtain  the  stresses  in  twisted  circular  cylinders  com¬ 
posed  of  the  same  materials. 

A  solution  is  also  found  for  the  stresses  in  a  thick  spherical 
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shell  expanded  by  an  internal  pressure  and  composed  of  a  neo-Hookean 
el astic-perfectly  plastic  material  obeying  any  isotropic  yield  condition. 
The  residual  stresses  resulting  from  unloading  from  an  elastic-plastic 
state  are  also  found. 
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CHAPTER  I 

PERFECTLY  PLASTIC  AND  ELASTIC-PERFECTLY  PLASTIC  CQNTINUA 

1 . 1  Introducti on 

A  branch  of  applied  mathematics  called  continuum  mechanics 
has  been  developed  which  attempts  to  describe  the  macroscopic  behavior 
of  deforming  materials  by  considering  idealized  mathematical  models. 

It  is  not  surprising  in  view  of  the  natural  division  of  materials  into 
solids  and  fluids,  that  the  two  most  well  established  branches  of  this 
field  are  classical  elasticity  and  hydrodynamics .  The  need  however, 
for  more  elaborate  theories,  valid  for  the  wide  range  of  applications 
of  materials,  and  to  describe  the  numerous  new  materials  which  have 
been  developed  has  accelerated  the  growth  of  continuum  mechanics  and 
has  led  to  several  new  branches. 

One  such  development  is  the  mathematical  theory  of  plasticity 
The  plastic  deformation  of  metals  is  the  basis  of  all  metal  forming 
processes  and  consequently  it  is  of  considerable  technological  interest 

It  is  observed  in  a  simple  tension  test  of  a  material  such  as 
steel  that  elastic  behavior  occurs  under  increasing  load  only  up  to  a 
certain  point  called  the  yield  point  after  which  plastic  flow  occurs 
with  continuing  deformation  resulting  with  little  increase  in  load. 

In  contrast  to  the  viscous  flow  of  fluids,  the  resistance  to  shear  de¬ 
formation  of  plastically  flowing  materials  appears  to  be  almost  inde- 
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pendent  of  the  rate  at  which  the  deformation  occurs. 

Early  work  in  plasticity,  such  as  that  of  Coulomb  [1]*  in 
1773,  dealt  mainly  with  the  failure  of  soils.  In  1864  Tresca  [2] 
published  the  results  of  some  experiments  on  the  extrusion  of  metals 
and  concluded  from  these  experiments  that  the  yield  point  is  reached 
when  the  maximum  shearing  stress  reaches  a  critical  value.  Soon  after, 
St.  Venant  [3]  suggested  that  for  plane  plastic  strain  the  principal 
axes  of  stress  and  plastic  strain  increment  coincide.  L£vy  [4]  in  1870 
extended  this  idea  to  three  dimensions  and  proposed  a  coaxial  relation 
between  the  stress  deviator  tensor  and  the  strain  increment  tensor. 

The  same  relation  was  suggested  independently  by  von  Mises  [5]  in  1913. 

In  1928  von  Mises  [6]  extended  his  theory  to  consider  a  perfectly  plastic 
solid  with  an  arbitrary  regular  yield  function  f(a.  .). 

*  sJ 

Hencky  [7]  proposed  a  so  called  total  deformation  theory  in 
1924  in  which  the  stress  and  the  total  plastic  strain  are  related 
rather  than  the  stress  and  the  plastic  deformation  rate  or  strain  in¬ 
crement  as  in  the  incremental  or  flow  theory  of  St.  Venant  and  L£vy. 

Valid  objections  have  been  raised  against  total  deformation  theories  by 
Handelman,  Lin,  and  Prager  [8]  and  by  Hill  [9].  These  theories  have 
often  been  used  however,  especially  in  Russia,  and  under  some  conditions 
they  give  satisfactory  results. 

The  incorporation  of  elastic  strains  in  the  incremental  or 
flow  theory  was  first  done  by  Prandtl  [10]in  1924  for  plane  problems  and 

*The  numbers  in  square  brackets  refer  to  references  listed  in  the 
Bibliography  at  the  end  of  this  thesis. 
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by  Reuss  [11]  in  1930  for  three  dimensional  problems. 

In  the  period  following  World  War  II  until  the  present, 
much  work  has  been  done  to  put  the  theory  of  plasticity  on  a  sound 
mathematical  and  physical  foundation.  Several  reference  books  have 
been  written  [12],  [13],  [14],  [15],  [16],  [17],  and  a  number  of  ex¬ 
tensive  review  papers  have  appeared  [18],  [19], 

Some  recent  developments  in  the  theory  of  elastic-plastic 
conti nua  have  involved  the  consideration  of  finite  deformations.  For 
example  the  kinematically  correct  form  of  the  Prandtl -Reuss  equations, 
valid  for  large  total  deformations  with  small  elastic  strains  was  first 
presented  by  Thomas  [15],  [20]  in  1955. 

Classical  elastic-plastic  theory,  including  the  Prandtl - 
Reuss  equations,  is  not  applicable  to  the  consideration  of  elastic- 
plastic  flow  in  the  presence  of  finite  elastic  strain.  One  example 
of  this  is  the  elastic-plastic  flow  of  a  metal  which  is  subjected  to 
an  hydrostatic  pressure  which  is  sufficient  to  produce  finite  elastic 
volume  change.  Such  pressures  are  produced  in  certain  explosive  forming 
processes . 

The  results  of  some  preliminary  experiments  at  the  Department 
of  Mechanical  Engineering,  University  of  Alberta,  indicate  that  several 
types  of  a  urethane  elastomer  (duPont  Adiprene)  exhibit  permanent  set 
after  finite  elastic  shear  strains.  Cubes  of  this  material,  of  5/8 
inch  side,  were  compressed  slowly  between  lubricated  parallel  flat  dies. 
Tests  conducted  in  the  temperature  range  70  -  80°F  showed  significant 
rate  effects  with  the  recovery  of  a  major  portion  of  the  residual 


3a 

deformation  occurring  over  a  period  of  several  days.  At  temperatures 
around  32° F  however  the  residual  deformation  which  occurred  after  red¬ 
uctions  of  approximately  30%  appeared  to  be  permanent.  The  elastic 
behavior  of  these  materials  is  approximately  neo-Hookean  and  volume 
measurements  of  deformed  specimens  showed  that  they  are  very  nearly 
incompressible . 

These  experiments  were  not  extensive  and  are  not  part  of  the 
work  of  this  thesis.  They  are  mentioned  however  because  they  indicate 
a  practical  and  important  area  for  more  extensive  experimental  work  which 
may  or  may  not  justify  some  of  the  assumptions  made  in  this  thesis. 

High  energy  rate  metal  forming  problems  and  the  possibility 
that  an  elastic-plastic  theory  may,  although  not  yet  experimentally 
verified,  provide  an  approximation  to  the  behavior  of  some  elastomers  at 
certain  temperatures  has  motivated  the  development  of  a  number  of  theories 
for  elastic-plastic  continua  which  may  undergo  finite  elastic  strains. 

One  of  the  first  works  which  may  be  applicable  to  such  problems 
appears  to  have  been  by  Sedov  [21].  A  theory  for  rate  independent  mater¬ 
ials  in  which  the  stress  is  a  functional  of  the  entire  deformation  hist¬ 
ory  has  been  presented  by  Pipkin  and  Rivlin  [22].  Isothermal  deformations 
only  are  considered.  A  theory  which  is  not  restricted  in  this  manner 
has  been  proposed  by  Green  and  Naghdi  [23],  [24].  It  is  modelled  on 
the  classical  flow  theory  and  makes  considerable  use  of  thermodynamic 
considerations.  Lee  [25],  [26],  [27]  has  also  suggested  a  flow  theory 
for  elastic-plastic  materials  under  finite  elastic  and  plastic  strains 
which  differs  from  that  of  Green  and  Naghdi.  The  most  recent  work  is 
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that  of  Freund  [28].  It  is  an  extension  of  the  work  of  Lee  and  introduces 
a  plastic  rate  of  deformation  tensor  which  differs  from  that  suggested 
by  Lee  and  that  developed  in  Chapter  II. 

In  this  thesis  the  kinematics  for  a  general  theory  for  materials 
which  may  undergo  finite  elastic  and  plastic  strains  is  discussed  in 
Chapter  II.  Following  the  discussion  in  section  3.7  the  development  of 
the  theory  is  restricted  to  materials  which  are  elastically  isotropic. 

An  argument  based  on  the  assumption  that  the  materials  considered  are 
elastically  incompressible  and  obey  Drucker's  postulate  [30]  leads  to  a 
general  plastic  flow  rule  from  which  two  specific  flow  rules  are  found, 
one  of  which  is  associated  specifically  with  materials  whose  elastic 
behavior  is  neo-Hookean.  A  simplification  which  results  if  plastic  iso¬ 
tropy  is  assumed  is  indicated  in  section  4.4.  Three  boundary  value  problems 
are  solved.  These  are  the  isothermal  simple  shear  of  a  cuboid,  torsion 
of  a  circular  cylinder,  and  expansion  of  a  spherical  shell,  the  respective 
deformations  being  homogeneous  and  rotational,  inhomogeneous  and  rotational, 
and  inhomogeneous  and  i rrotational .  The  materials  considered  in  the  prob¬ 
lems  are  isotropic,  incompressible,  and  exhibit  neo-Hookean  elastic  behavior. 


1 .2  Constitutive  Equations  for  an  Hookean  Elastic  Material 

Let  u.  be  the  components  of  the  displacement  vector  of  a 
particle  of  an  Hookean  elastic  material,  referred  to  a  fixed  Cartesian 
coordinate  system.  The  small  strain  tensor  referred  to  this  coordinate 
system  has  components 


i  3u.  3u . 

e  =  -  (— 1  +  — J-) 

6i j  2  l3Xi  3x.;  ' 


(1.2.1) 


The  constitutive  equation  for  an  Hookean  elastic  material  is 
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°1j  '  l+v  ei j  +  (l+v)(l-2v)  ekk  6ij  (1-2‘ 

where  are  the  components  of  the  Cauchy  stress  tensor,  E  is  Young's 
modulus,  and  v  is  Poisson's  ratio  [29]. 

Substituting  the  components  of  the  stress  deviator 


3  °kk 


and  the  strain  deviator 


3  ekk 


in  equation  (1.2.2)  gives 


e*  =-L 
ij  2p 


°ij  •  ekk 


=  -  P/K 


(1 ,2.3a  ,b) 


F  £ 

where  y  =  ^ ( -[  +v )  is  the  modulus  of  rigidity,  K  =  3( y_2v)  1>s 

modulus  and  p  =  For  an  incompressible  material  el.  =  e.  .  and 

only  the  first  of  equations  (1.2.3)  is  required.  The  hydrostatic  part 


of  the  stress  is  then  determined  from  the  equations  of  motion. 


1 .3  Constitutive  Equations  for  a  Rigid  Perfectly  Plastic  Material 

The  rigid  perfectly  plastic  material  is  a  hypothetical  material 
which  remains  rigid  under  loading  until  the  stresses  satisfy  a  yield 


' 


6 


condition.  Then  irreversible,  rate  independent  deformation  occurs  at 
constant  stress.  The  model  is  a  reasonable  approximation  to  some 
structural  materials  such  as  steel  which  have  a  well  defined  yield 
point  and  for  which  the  modulus  of  rigidity  is  large  compared  to  the 
yield  stress.  It  is  not  suitable  however  for  problems  in  which  the 
plastic  strains  are  constrained  to  be  of  the  same  order  of  magnitude  as 
the  elastic  strains. 

Fundamental  to  the  theory  is  the  concept  of  a  yield  condition. 
It  is  postulated  that  there  exists  a  function  of  the  stress  and  temper¬ 
ature,  f(a..,T),  called  the  yield  function  which  has  the  properties 
'  J 

f  <  0  , 


dj*?)  =  0  if  f  <  0  or  f  =  0  and  df  <  0  , 

'  J 


and 


d!1?)  |  >  0  if  f  =  0  and  df  =  0  , 
’  j 


where  d!^  are  the  components  of  the  plastic  rate  of  deformation  tensor 

*  vJ 

given  by 


/  -i  3v.  3v. 

d(P)  =  1  ( L  +  i) 

°i j  2  '3x.  3x.'  * 

J  1 

v^  being  the  components  of  velocity  at  a  point. 

In  nine  dimensional  stress  space  the  yield  condition  f(a.  .,T)  =  0 

^  J 

defines  a  closed  surface.  Deformation  is  possible  only  under  stress  systems 


. 


■ 


which  lie  on  the  surface,  and  stress  points  outside  the  surface  are  not 
admissible. 

Consider  a  body  in  equilibrium  under  an  homogeneous  stress 
system  ct?j  which  is  inside  or  on  the  yield  surface.  Suppose  that  an 
independent  external  agency  applies  and  removes  isothermal ly  an  addi¬ 
tional  stress  to  the  body.  Drucker's  postulate  [30]  states  that  the 
work  done  by  the  external  agency  during  the  complete  cycle  is  non¬ 
negative  for  a  perfectly  plastic  material.  That  is,  the  external 
agency  cannot  extract  useful  work  from  the  system  by  applying  and  re¬ 
moving  additional  stresses.  Drucker's  postulate  defines  a  class  of 
material  models  which  have  been  found  to  provide  an  acceptable  ap¬ 
proximation  for  many  structural  materials. 

Two  important  results  follow  from  this  postulate.  Let  the 

stress  point  a?,  lie  inside  the  yield  surface.  With  the  addition  and 

^  J 

removal  of  the  external  agency  the  stress  point  follows  a  closed  path 
in  stress  space.  Suppose  that  during  the  time  interval  [t,  t  +  6t] 
the  stress  path  lies  on  the  yield  surface.  That  is,  during  this  time 
interval,  the  body  has  yielded.  The  work  per  unit  volume  done  by  the 
total  stress  system  during  the  complete  cycle  is 


t+6t 


W  = 


t 


dt  , 


since  for  a  rigid  plastic  material  deformation  can  occur  only  when  the 
stress  point  lies  on  the  yield  surface  and  in  general  for  an  elastic- 


’  rtteq  baaofa  s  awoHtrt  Jnroq  aestJi  JiiJ  1CM««  remsSxs  srtt  4o  svoma-i 

" 
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plastic  material  the  net  work  done  on  the  elastic  strains  in  a  closed 
cycle  is  zero. 

Similarly  the  work  done  per  unit  volume  by  the  original 
stress  a?,  during  the  complete  cycle  is 

*  J 


0T  = 


t+6t 

°ijdlj)dt 


Drucker's  postulate  states  that 


t+6t 


W  -  or  = 


(oij-°ij)dij)dt-° 


It  follows  that 


(a .  .-a?  .  )d >  0  , 
ij  iJ  iJ  - 


(1.3.1) 


where  a?,  is  any  stress  point  lying  on  or  inside  the  yield  surface  and 

J 

a.  •  is  any  stress  point  on  the  yield  surface  with  an  associated  plastic 

'  J 

(P) 

rate  of  deformation  tensor  d..'  . 

i  J 

Equation  (1.3.1)  which  is  known  as  the  maximum  work  inequality 

implies  that  the  yield  surface  is  convex  and  that  a  vector  in  stress 

(P) 

space  with  components  d..'  ,  multiplied  by  a  suitable  scalar  to  give 

■  3 

the  units  of  stress,  is  normal  to  the  yield  surface  at  the  point  ck  . 
or  lies  between  adjacent  normals  if  ck  .  is  at  a  corner  of  an  irregular 
yield  surface. 


'  ■  *  •  ■  '  "  ■  i 
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Thus  for  a  regular  yield  function 


(1.3.2) 


where  X  is  a  non-negative  scalar  factor  of  proportionality  and  A  =  0 


if  f  <  0  or  f  <  0  .  The  flow  rule  (1.3.2)  implies  that  if  yielding  is 


not  influenced  by  the  hydrostatic  part  of  the  stress  tensor  the  plastic 


deformation  is  isochoric,  that  is  dfP^  =  0. 


For  an  isotropic  rigid  plastic  material  the  form  of  the  yield 


function  at  a  point  must  be  independent  of  the  orientation  of  the 
material  element  considered  so  that  the  yield  function  depends  only  on 
the  invariants  of  the  stress  tensor.  Furthermore  the  experimentally 
observed  fact  that  yielding  is  almost  independent  of  moderate  hydro¬ 
static  pressure  [31]  results  in  the  usual  assumption  that  for  an  iso¬ 
tropic  solid  f  is  a  function  of  the  stress  deviator  invariants 


and 


In  general,  the  flow  rule  associated  with  the  yield  function 
f  for  an  incompressible  isotropic  perfectly  plastic  material  is 


bOK  ortJ  ’to  imo^  M#b  orqo-  Joe  ne  W 


■ 


10 


For  the  von  Mises  yield  condition 

1  °i j°i j  "  k2  =  0  >  (1.3.4) 

where  k  is  the  yield  stress  in  pure  shear,  the  flow  rule  (1.3.3)  be¬ 
comes 


3 


from  which 


d!P)d<P)=  X2a'..ai. 
1J  1J  1J  1J 


(1.3.5) 


From  equations  (1.3.4)  and  (1.3.5)  the  following  expression  for  A  is 
obtained , 


■d!p>d<p) 

1J  1J 


and  thus 


,(p)_ 


rs 


A 


(P)d(P? 


II 


As  mentioned  by  Thomas  [32],  there  is  no  one-one  relationship 
between  the  stress  deviator  and  the  plastic  rate  of  deformation  tensors 
unlike  the  constitutive  equations  for  a  viscous  fluid. 


n 


1 .4  Prandtl-Reuss  Equations  for  Small  Total  Deformation 

The  classical  elastic-plastic  material  is  one  for  which  the 
total  deformation  rate  may  be  written  as  the  sum  of  a  plastic  and  an 
elastic  part  which  are  obtained  respectively  from  the  rigid  plastic 
flow  rule  and  the  material  derivative  of  the  Hookean  elastic  stress 
strain  law. 

Suppose  that  during  a  time  increment  dt  the  increment  of 

stress  at  a  point  in  an  elastic-plastic  material  is  da.,  referred  to 

^  J 

a  fixed  Cartesian  coordinate  system  ox.  and  denote  by  du.  the  increment 
of  displacement  which  occurs  during  this  time  increment.* 

The  total  strain  increment  then  is 


1  3(du.)  3(du.) 

de  =  -  ( _ 1—  +  _ 1-1 

Qei j  2  1  3x .  3x .  ' 

J 


(1.4.1) 


It  is  assumed  that  de. .  may  be  written  as  the  sum  of  an  elastic  and  a 

^  J 

plastic  strain  increment**  so  that 


de 


(e) 

'ij 


de(P> 

IJ 


and  thus 


d..  .  d<e> 


TJ 


1  J 


d(P} 

IJ 


where  d. 


(e) 


ij 


and  d(P^ 

*  J 


are  the  elastic  and  plastic  rate  of  deformation 


*  The  components  du..  are  given  by  du^  =  v^ dt  where  v.  is  the  velocity 
at  the  point  considered. 

**A1 though  de.  .  must  be  derivable  from  an  incremental  displacement  field 
iJ  (e)  (d) 

this  is  not  necessarily  true  for  de:  .  and  de_.^'  . 

I  J  •  J 
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tensors  respectively. 

An  expression  for  the  elastic  rate  of  deformation  tensor  is 
obtained  by  taking  the  material  derivative  of  equations  (1.2.3)  to  give 


=  J-  a! .  -  Jr  p  6 .  . 
l  J  2  y  ij  3K  K  i  j 


The  plastic  rate  of  deformation  tensor  is  given  by  equation 
(1.3.2)  and  thus  the  total  deformation  rate  tensor  during  loading  is 


j  -j:,  j:  *  I  •»  3f 

di j  '  2y  °ij  "  3K  p  5ij  X 


30  ij  ' 


In  particular  for  an  elastic-perfectly  plastic  material  obeying  the 
von  Mises  yield  condition  (1.3.4)  the  constitutive  equation  is 


_1_ 

2p 


±_ 

3K 


(1.4.2) 


An  expression  for  X  may  be  obtained  by  multiplying  equation 


(1.4.2)  by  and  contracting.  This  gives 


d .  .o  \  . 
X  =  _LL_L1 

A  2 

2k 


with  X  =  0  if  a!,  a!  .  <  2k^  or  if  d..  a!.  <  0  . 

•  J  ^  J  1  J  •  J 

Equations  (1.4.2)  are  called  the  Prandtl-Reuss  equations 
and  for  an  incompressible  material  they  become 


1  Lpill 


13 


2y 


a .  . 

ij 


Aaij 


(1.4.3) 


1 .5  Invariance  Requirements  for  Constitutive  Equations 

The  constitutive  equations  of  a  material  are  essentially  a 
mathematical  description  of  the  deformation  properties  of  the  material 
and  must  satisfy  certain  invariance  requirements  to  be  physically 
meaningful . 

The  first  requirement  is  called  coordinate  invariance.  The 
deformation  characteristics  of  a  material  must  be  independent  of  the 
coordinate  system  used  to  describe  the  system.  This  requirement 
is  satisfied  by  representing  correctly  all  physical  quantities  by 
their  components  as  scalars,  vectors,  and  tensors  referred  to  some 
coordinate  system. 

The  second  invariance  principle  may  be  stated  in  two  slightly 
different  but  equivalent  forms  called  the  principle  of  material  in¬ 
difference  or  material  objectivity  and  the  principle  of  isotropy  of 
space  [33],  [34],  [35]. 

Consider  a  body  whose  motion  is  given  by  (Xw,t)  referred 
to  a  coordinate  system  in  the  reference  frame  S.  The  quantities  x.  are 
the  coordinates  at  time  t  of  a  material  particle  which  occupied  the  point 
with  coordinates  at  some  reference  time  t  .  Consider  the  following 
equations , 


xi(XK,t)  =  C.(t)  +  Q i j ( t )  Xj(XK,t)  , 


(1.5.1) 


. 

. 


. 
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where  Q-k  Qjk  =  6.^.  for  all  time  t. 

Two  interpretations  may  be  given  to  this  equation  [33]. 

Firstly  it  may  be  thought  to  relate  the  same  motion  of  a  body  as  seen 

by  observers  in  two  different  reference  frames,  S  and  S,  moving 

relative  to  one  another.  For  this  interpretation  the  transformation 

Q..(t)  belongs  to  the  full  orthogonal  group,  that  is  transformations 
^  J 

for  which  det[Q.  .]  =  -  1  are  included. 

*  J 

Equation  (1.5.1)  may  also  be  interpreted  as  relating  two 
motions  of  a  body,  differing  only  by  a  rigid  body  motion,  as  seen 
by  a  single  observer.  To  be  meaningful,  Q. .  must  then  belong  to  the 

•  J 

group  of  proper  orthogonal  transformations,  that  is  det[Q..]  =  +1. 

*  J 

Returning  to  the  first  interpretation  of  equation  (1.5.1), 
the  following  definition  is  useful.  Quantities  whose  transformation 
from  reference  frame  S  to  S  depends  only  on  the  relative  orientation 
of  the  two  reference  frames  and  not  on  any  other  aspect  of  their  re¬ 
lative  motion  are  said  to  be  frame- indifferent  or  objective.  Thus 
frame-indifferent  scalars,  vectors,  and  tensors  transform  as 

A  =  A 


vi  =  Qij(t) 


kl  * 


and 


bi j  =  Qik^  Qji*1)  b 


. 
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respectively  for  all  time  t. 

Since  the  response  of  a  material  to  some  loading  system 
must  appear  the  same  to  all  observers  regardless  of  their  motion,  con¬ 
stitutive  equations  must  be  form  invariant  under  a  change  of  reference 
frame.  Thus  in  general  if  the  constitutive  equation  for  a  material 
as  observed  from  reference  frame  S  is  of  the  form 

ai j  =  Hij^xk^XK,t^  ’ 

where  the  tensor  components  H.  .  are  functionals  of  the  entire  defor- 

*  J 

mation  history  of  the  body,  then  the  constitutive  equation  for  the 

<v 

same  material  as  observed  from  the  reference  frame  S  must  be 

aij  =  Hij^xk^XK#t^  ’ 

where  a...  =  Qi(<(t)  Qj-j(t)  and  xk(XK,t)  is  given  by  equation  (1.5.1). 
This  same  result  follows  if  the  second  interpretation  of  equation 
(1.5.1)  is  used. 

It  will  now  be  shown  that  the  Prandtl-Reuss  equations  given 
in  the  previous  section  do  not  satisfy  the  principle  of  material  in¬ 
difference.  The  results  one  obtains  using  the  Prandtl-Reuss  equations 
for  any  particular  problem  must  be  considered  in  view  of  this  since 
the  use  of  these  equations  in  situations  where  rotations  are  significant 
may  result  in  unacceptable  errors. 

Taking  the  total  time  derivative  of  equation  (1.5.1)  gives 


. 


' 
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*>* 


c.(t)  +  Q-j  j(  t) 


Defining  the  antisymmetric  tensor 


n. . 

ij 


by 


j  ~  Q j k  Q j |<  ”  “  j  (1.5.2a,b) 

the  velocity  gradient  tensor  observed  from  reference  frame  S 
is 


3v.  9v 

-  Q  ^  l  Q  ^ i  +  ^ .  (1.5.3) 

ax.  1K  Jl  ax,  1J 

J 


Forming  the  symmetric  and  antisymmetric  parts  of  the  velocity  gradient 
tensor  gives 


and 


!ik 


qji 


wkl 


+  Q 


1J 


* 


which  are  the  stretching  or  rate  of  deformation  and  the  spin  or 
vorticity  tensors  respectively.  Thus  the  rate  of  deformation  tensor 
is  frame- indifferent  and  the  spin  tensor  is  not. 

Consider  the  total  time  derivative  of  the  stress  tensor 


°ij  Qik  Qjl  akl  ' 


- 
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This  leads  to 


:ij  =  QikQjlakl 


QikQjlakl 


QikQjlakl 


(1.5.4) 


which  shows  that  the  material  derivative  of  the  stress  tensor,  and 
hence  the  stress  deviator  tensor,  is  not  frame-indifferent  and  thus 
the  Prandtl -Reuss  equations  (1.4.3)  are  not  form  invariant  under  a 
change  in  reference  frame.  If  the  material  derivative  of  stress  in 
equation  (1.4.2)  was  replaced  by  a  frame-indifferent  stress  rate  then 
it  follows,  since  d. .  and  a  I.  are  frame- indifferent,  that  equation 

T  J  ^  J 

(1.4.2)  would  be  form  invariant  under  changes  of  reference  frame. 

From  equations  (1.5.2)  and  (1.5.3)  it  is  found  that 


3v.  9v 

Q  .  -  Q. 

3xr  4rk  wir  3xk 


3v. 

3x.j  ^rk  +  r  3xr  ' 


(1.5.5) 

(1.5.6) 


Substituting  equations  (1.5.5)  and  (1.5.6)  into  equation 
(1.5.4)  gives 


a 


o 


9v. 
rj  Sx" 


9v. 


"  Q,°  bQ-n  ( 


3  V. 


9v 


ir  3x 


-a. 


i kyjl v  kl  wrl  9x 


-  a 


kr  3x. 


)  ,  (1.5.7) 


4  +„ 

a .  .  +  a  .  7T — 

ij  rj  9x 


9vr  .  9vr 


9v 


+  5ir  3^  =  QikQjl(okl+arl  3^+  akr  3^  >  <1-5-8) 


and 


respectively.  Addition  of  equations  (1.5.7)  and  (1.5.8)  and  use  of 
the  definition  of  the  spin  tensor  gives 


a.  .  +  co  .a  .  +  w  .a.  =  Q  °  ■  Q  •  i  (a,  t+go  ,  a  n  +co  na,  )  .  (1.5. 

ij  r i  rj  rj  ir  yikyjlv  kl  rk  rl  rl  kr  v 

Using  the  notation 


6a 


6t 


U  = 


9v 


l 


3vi 


aij  "  arj  9xr  "  air  9xf  * 


9v  8v 

i/  r  r 

7T7-  a.  .  =  a.  •  +  a  .  t —  +  a.  -r —  , 
ft  ij  ij  rj  9x.  lr  9x . 


and 


8.  .  -  a.  .  +  a  .w  .  +  a.  a)  .  , 

i  j  i  j  rj  n  l  r  r  j  * 


equations  (1.5.7),  (1.5.8),  and  (1.5.9)  become 


6a. 


■fit*1  "  kQjl  6t  °kl  ’ 


V  ~  _  n  a  V 

Vt  aij  yik4jl  Vt  ak1  ’ 


and 


o 


a .  .  =  Q  . ,  Q  .  ,a,  . 

ij  yikxjl  kl 


This  indicates  that  these  stress  rates,  the  Oldroyd  [36],  the  Cotter- 
Rivlin  [37],  and  the  Jaumann  [38]  stress  rates  respectively,  are  all 


■ 
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frame-indifferent  and  may  therefore  be  used  in  constitutive  equations. 

Other  stress  rates  may  be  formed  from  any  of  the  above 
stress  rates  by  addition  of  terms  which  are  frame- indifferent. 

1 .6  The  Kinematically  Correct  Form  of  the  Prandtl-Reuss  Equations 

As  shown  in  the  previous  section  the  use  of  the  material 
derivative  of  the  stress  or  stress  deviator  tensor  in  constitutive 
equations  violates  the  principle  of  material  indifference.  Thus  one 
of  the  frame-indifferent  stress  rates  must  be  used  instead. 

It  has  been  suggested  by  Prager  [39]  that  the  most  suitable 
definition  of  stress  rate  for  the  theory  of  plasticity  is  the  Jaumann 
derivative,  since  only  it  has  the  property  that  its  vanishing  implies 
that  the  stress  invariants  are  stationary.  Since  the  yield  condition 
for  an  isotropic  material  depends  on  the  invariants  of  the  stress  de¬ 
viator  and  since  any  observer  who  observes  zero  rate  of  change  of 
stress  should  observe  no  change  in  the  value  of  the  yield  function, 
the  choice  of  the  Jaumann  stress  rate  seems  justified. 

The  following  physical  interpretation  may  be  given  to  the 
Jaumann  derivative  which  adds  to  its  desirability  for  use  in  constitutive 
equations.  The  Jaumann  derivative  of  the  stress  tensor  at  a  point 
in  a  deforming  medium  is  the  material  derivative  of  the  stress  referred 
to  a  Cartesian  coordinate  system  which  instantaneously  coincides  with 
the  fixed  coordinate  system  but  which  participates  in  the  rotation 
of  the  material  particle  at  that  point.  Such  a  coordinate  system  has 
been  called  a  kinematically  preferred  coordinate  system  by  Thomas  [40] 


■ 
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and  the  Jaumann  derivative  is  sometimes  referred  to  as  the  co-rotational 
deri vati ve. 

The  kinematically  correct  or  frame-indifferent  form  of  the 
Prandtl-Reuss  equations  using  the  Jaumann  stress  rate  were  first 
developed  by  Thomas  [20]. 

After  substitution  of  the  Jaumann  stress  rate  for  the  ma¬ 
terial  derivative  of  stress,  the  Prandtl-Reuss  equations  become 


2u 


ij 


3K 


P  6 


1J 


Xal. 


9 


since 


p  -  p  • 


In  particular  for  an  incompressible  material 


JL 

2p 


°ij  +  Xaij 


o 


(1.6.1) 


This  form  of  the  Prandtl-Reuss  equations  is  valid  for  large 
total  deformations  since  the  use  of  frame-indifferent  tensors  allows 
for  finite  rotations  of  the  material  elements.  They  are  valid  how¬ 
ever,  only  for  small  elastic  strains,  that  is  elastic  strains  which 
may  be  described  by  the  small  strain  tensor  (1.2.1). 

The  constitutive  equations  (1.6.1)  have  been  used  by  Haddow 
and  Danyluk  [41],  [42],  [43]  to  solve  the  problems  of  the  plane  and 
axi symmetric  flow  of  an  incompressible  elastic-perfectly  plastic 
material  in  a  converging  channel. 


' 
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CHAPTER  II 

KINEMATICS  FOR  A  LARGE  STRAIN  ELASTIC-PLASTIC  THEORY 

2 . 1  Introduction 

In  this  chapter  the  kinematics  of  a  large  strain  elastic- 
plastic  theory  are  discussed  and  frame- indifferent  tensor  quantities 
are  developed  which  may  be  used  in  constitutive  equations  without 
violating  the  principle  of  material  indifference.  The  only  assumption 
made  is  that  both  finite  elastic  and  plastic  strains  may  be  possible 
in  the  materials  considered. 

Three  configurations  of  a  body  composed  of  an  elastic-plastic 
material  are  defined.  These  are  the  initial  unstrained  and  unstressed 
configuration  of  the  body  at  a  uniform  reference  temperature,  the 
unstressed  plastically  strained  configuration  of  the  body  at  the 
reference  temperature,  and  the  stressed  configuration  of  the  body  in 
which  elastic  and  possibly  plastic  strains  may  be  present  due  to  ex¬ 
ternal  loads  and  body  forces.  These  configurations  of  the  body  will 
be  referred  to  as  C.l,  C.2,  and  C.3  respectively.  If  no  plastic  strains 
have  occurred  then  C.l  and  C.2  differ  at  most  by  a  rigid  body  motion. 

The  use  of  three  configurations  to  consider  finite  elastic-plastic 
problems  has  been  adopted  by  Sedov  [21],  Backman  [44],  and  recently 
by  Lee  [27]. 

If  the  deformation  from  C.l  to  C.3  is  non-homogeneous  any 


;  fMfiw 

DU26fq  on  .y.rovttosq*!^  £.3  btis  ,S  3  .  ■  86  °:  -  9d 
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plastic  strains  which  result  will  in  general  not  satisfy  the  compatibility 
conditions  for  finite  strain  [45],  Thus  residual  stresses  will  result 
from  removal  of  the  external  loads  and  body  forces  required  to  maintain 
the  body  in  C.3  and  the  space  corresponding  to  C„2  will  then  be  non- 
Euclidean  and  C.2  will  not  be  a  configuration  which  the  body  may  occupy 
physically.  The  configuration  C.2  is  then  thought  of  as  the  configuration 
of  the  body  which  results  from  isolating  each  material  element  from  the 
surrounding  material  and  unloading  it  elastically  and  then  returning 
it  to  the  reference  temperature. 

A  fourth  configuration  of  the  body  C.4,  which  is  used  in 
Chapter  VII  results  from  removing  all  external  loads  and  body  forces 
from  C.3  and  returning  the  body  to  the  uniform  reference  temperature. 

If  the  plastic  strains  satisfy  the  compatibility  conditions  then  C.2 
and  C.4  coincide,  otherwise  residual  stresses  are  present  in  C.4. 


2.2  Deformation  Gradients  and  Strain  Tensors 

Let  the  coordinates  of  a  particle  in  C.l,  C.2,  and  C.3,  re¬ 
ferred  to  a  fixed  Cartesian  coordinate  system,  be  X^,  xa»  and  x.  re¬ 
spectively.  Henceforth  lower  case  Latin  subscripts  refer  to  C.3, 
lower  case  Greek  subscripts  to  C.2,  and  upper  case  Latin  subscripts 
to  C.l. 

The  motion  of  the  body  is  given  by 


xk(XK,t)  , 


(2.2.1) 


. 

' 
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or 


xk  =  xk^a,t^  and  *a  =  ^(XK*t)  > 


(2.2.2a, b) 


where  in  general  the  functions  x^fX^.t)  and  x^(x  »t)  have  different 
forms  o 

Continuity  requires  that  the  mapping  (2.2.1)  from  C.l  to  C.3 
have  a  unique  inverse  in  the  neighbourhood  of  any  given  particle.  It 
follows  from  the  implicit  function  theorem  [46]  that  x^X^t)  is  con¬ 
tinuous  and  possesses  first  order  partial  derivatives  with  respect  to 
X^  in  the  neighbourhood  of  a  particle  and  that  det  (ax^/SX^)  does  not 
vanish  in  that  neighbourhood.  These  conditions  are  not  necessarily 
true  for  the  mappings  from  C.l  to  C.2  and  C.2  to  C.3  since  C.2  may 
not  be  a  physically  attainable  configuration  of  the  body. 

Let  dS ,  di ,  and  ds  be  an  infinitesimal  element  of  length 
in  C.l,  C.2,  and  C.3  with  components  dX^,  dx ,  and  dx^  respectively, 
referred  to  a  fixed  Cartesian  coordinate  system.  Then 


dx.  s  F.^dX,,  -  F.  dy 
i  iK  K  ia  *01 


(2,2,3a,b) 


and 


dx  =  F(P'dX,, 

Aa  aK  K 


(2.2.4) 


3x 


( 6  ) 

where  F.^  =  is  the  total  deformation  gradient  and  F^y  and  F 


K 


(p) 

aK 


are  the  elastic  and  plastic  deformation  gradient  tensors  respectively 


These  tensors  are  related  by 


. 


. 


V 
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If  the  space  corresponding  to  C,2  is  Euclidean  then 


and 


F(p>  =  lOs. 


aK 


■a 


If  this  space  is  non-Eucl idean  the  above  partial  derivatives  do  not 
exist  but  equations  (2,2*3)  and  (2*2,4)  are  still  valid  with  the 


matrices  F?^  and  representing  local  1  inear  mappings .  Although 


the  remainder  of  the  discussion  in  this  chapter  is  restricted  to  con¬ 
figurations  C.2  which  are  Euclidean,  the  results  remain  valid  in 
general  provided  velocity  and  deformation  gradients  with  respect  to 
C,2  are  not  expressed  as  partial  derivatives. 

The  assumption  from  classical  plasticity  that  the  plastic 
volume  change  is  zero  is  retained  here  and  thus 


(2,2,5) 


Furthermore 


detf^1)  =  det^)  =  — 


3x»  3x.  p. 


9 


. 


25 


where  p-j  and  p^  are  the  densities  in  C.l  and  C.3  respectively. 
From 


and 

and  equations 
tained , 


dS2  =  dXKdX„  , 

d"2  =  %%  1 

ds2  =  dx«dx.  , 

(2.2,3a,b)  and  (2.2.4)  the  following  results  are  ob- 


and 


where 


and 


ds2  -  dS2  =  2E|<LdX|<dXL  , 


ds2  -  d<52  =  2E(®V  dXo> 

aB  a  3 


di2  -  dS2  =  2E*P-dXKdXL  , 


1  3x.  3x. 

EKL  =  2  (3X^  IX^  '  6KL-  ’ 


/v  ,  3x.  9Xo 

E(e)  =  1  ( _ ] _ L  _  6  ) 


a 3  2  3xa 


a3J 


(p)  -  1  x  \ 

"  °KL  9 


EKl/  "2  ^3Xk  3X. 


(2 


2.6a  ,b  ,c) 


are  the  total,  the  elastic,  and  the  plastic  Lagrangian  strain  tensors 


. 
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referred  to  C.l,  C.2,  and  C « 1  respectively . 


2.3  Strain  Rate  and  Rate  of  Deformation  Tensors 

The  material  derivative  of  equation  (2.2.6a),  that  is  the 
time  derivative  holding  constant,  is 


3x .  9x . 


EKL  ■  di j  3Xk  3Xl  ’ 

(2.3.1) 

where 

i  3V.  3v. 

d  =  1  (  1  +  — 1) 

ai j  2  ^3x»  9x  = ' 

J 

(2.3.2) 

is  the  total  deformation  rate  tensor  and  v^  =  x^  is  the  velocity  of 
a  particle  in  C.3. 

(P) 

Similarly  material  differentiation  of  E^y  gives 


where 


;(p)  -  „(p) 

lKL  a 3  3Xk  9Xl  • 


AP) 

a  3 


t  3v 
I  /  a 


3v 


-  ( 

2  ^X 


3 


3 


9X, 


a 


is  a  plastic  rate  of  deformation  tensor  and  v  =  x  is  the  velocity  of 

v  a  Aa  J 

a  particle  in  C.2.  Equation  (2.2c5)  which  follows  from  the  assumption 
of  zero  plastic  volume  change  implies  that 


aa 


=  0  . 


27 


The  material  derivative  of  the  equation  (2.2.6b)  is* 


. /  \  ax.  ax.  n  sv  ax.  ax.  i  av  ax.  ax. 

c(e)  -  a  _ I  _JL  _  1  _X _ 1 _ 1 . 1  _X _ 1 _ 1 


Eai'  =  dij  ajj  '  2  <  3X^  3Xg  ‘  2  3^  3x^  ‘  (2’3'3) 


The  velocity  gradient  dv  /a xft  which  appears  in  equation 

Oi  lo 

(2.3.3)  may  be  written  as 


where 


+  03 


(P) 

aB  * 


(P)  .  1  ,^a 
aB  2  l3Xg  3xa 


(2.3.4) 


is  the  plastic  spin  tensor.  Configuration  C.2  is  defined  only  to 

(P) 

within  a  rigid  body  rotation  and  the  plastic  spin  tensor  is  not 

c(e)  (P) 

uniquely  defined.  The  strain  rate  E^'  depends  on  ur^'  and  consequently 
it  is  an  unsatisfactory  definition  of  elastic  strain  rate.  In  fact  it 
can  have  non-zero  components  when  the  body  is  undergoing  a  rigid  body 
motion  with  neither  elastic  nor  plastic  deformation  occurring.  The 
necessary  and  sufficient  condition  for  a  rigid  body  motion  is 


=  0 


throughout  the  body  and  the  condition  for  no  plastic  deformation  is 


*In  taking  this  material  derivative  it  is  noted  that  the  coordinates 
Xa  are  functions  of  time. 


. 

- 
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throughout  the  body.  If  both  conditions  are  met  then  no  elastic  de- 

(P) 

formation  is  occurring  but  is  in  general  not  zero  for  a  rigid 

body  motion  of  C.2  and  it  can  be  deduced  from  equations  (2.3.3)  and 

(2.3.4)  that  is  non-zero  when  d^.  =  0  and  dj^  e  0  if  f  0 

f  Pi 

An  elastic  strain  rate  tensor  which  does  not  depend  on  is  ob¬ 


tained  as  follows 


The  plastic  deformation  gradient  may  be  written  as 


_  R(P)y(P) 

SXK  '  V  ULK  * 


,(P) 


(P) 


where  J  is  a  proper  orthogonal  tensor  and  is  a  symmetric 


aL 

positive  definite  tensor  such  that 


u(P)u(P)  - 

UKM  ULM  "  8Xk  3Xl  * 


The  deformation  from  C.l  to  C.2  consists  of  a  pure  stretch  deformation 

(Pi  (Pi 

given  by  U^'  followed  by  a  rigid  body  rotation  given  by  .  There- 

( s ) 

fore  unlike  Ev  the  tensor 

a3 


E 


(e) 

KL 


=  R 


(P)n(P)r(e) 

aK  3L  a3 


(2.3.5) 


is  independent  of  any  rigid  body  rotation  of  C.2. 

The  material  derivative  of  equation  (2.3.5)  is 


j  yd  rievfg 

' 

I  frgt-r  y/is  Jnsbnsq -tool’  zt 
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r(e)  =  R(P)R(P)p(e)  .  E(e)  E(e)0 
KL  aK  8L  aB  4 KM  ML  KM  4 ML 


(2.3.6) 


where 


q  =1  (u(P)u(P> 

^KL  9  'U|/M  Umi 


KM  ML 


_1  -  u(p)_1u(p)) 

UKM  UML  ' 


(2,3.7) 


is  an  antisymmetric  tensor  and 


E^  =  +  E(eh(P)  -  JPV® 

a3  a3  ay  y3  ay  y3 


(2.3.8) 


is  the  co-rotational  derivative  of  E that  is  the  material  derivative 

referred  to  a  Cartesian  coordinate  system  which  participates  in  the 

rotation  of  the  material  in  C.2  at  the  point  considered  and  which  in- 

^  ( e ) 

stantaneously  coincides  with  the  fixed  coordinate  system.  Thus  E^y 
is  independent  of  the  spin  of  C.2, 

Using  equations  (2,3.3)  and  (2,3,4)  the  co-rotational  deri- 

(  p  ) 

vative  of  Ev  '  becomes 

a3 


|(e)  ,  ;(e)  +  (P)  +  l^i^i  (P) 

aB  aB  2  d\  %  yS  2  ^  3Xg  Y“ 


ij  3xa  SXg  "  2  ya  ZXy  3Xg 


1  h 

2  dyB 


3x.  3x. 


(2.3.9) 


An  elastic  rate  of  deformation  tensor  which  is  independent 


of  the  spin  of  C.2  is  given  by 


. 

. 

2tJ- 


. 


1  v J I  ■ . S  HI  } 


30 


,(e)  _  ^a^f(e) 
ij  "  ax.  ax.  a3  9 

\J 


and  from  equation  (2.3.9) 


d(e>  -  d 

ij 


. .  -  d* . 
i  J  iJ 


(2.3.10) 


where 


d*  =  i(!^d(p)  !!i  +  !^d(p)  ^i) 

ij  2  vax.  a3  3X  3x.  a3  9x 
°  l  Aa  j  Aa 


(2.3.11) 


The  tensor  d* .  is  independent  of  the  rotation  of  C.2  and 
i  J 

the  spin  of  C.2  and  C.3.  It  has  the  property  that  dj^  =  0  for  elastic- 
plastic  flow  with  no  plastic  volume  change.  Thus  d*.  is  a  possible 

1  J 

definition  for  plastic  rate  of  deformation  and  it  is  shown  in  section 
3.7  that  it  is  a  suitable  definition  for  plastic  rate  of  deformation 
for  an  elastic-plastic  material  which  is  elastically  isotropic. 
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CHAPTER  III 

HYPERELASTICITY  AND  THERMOELASTIC  CONSIDERATIONS 

FOR  AN  ELASTIC-PLASTIC  MATERIAL 

3.1  Introduction 

The  constitutive  equations  (1.2.2)  for  the  classical  Hookean 
elastic  material  are  based  on  a  rather  restricted  definition  of 
elastic  behavior.  In  fact  Hookean  elastic  behavior  is  a  special 
case  of  a  much  broader  concept  of  elasticity;  the  existence  of  a 
stress  free  natural  state  and  the  unique  dependence  of  the  stress 
on  the  deformation  gradients  relative  to  the  natural  state  and  on 
temperature  An  elastic  material  possesses  a  perfect  memory  of  the 
natural  state.  In  this  chapter  elastic  materials  [47]  which  possess 
an  elastic  potential  or  strain  energy  function  are  considered  and 
these  materials  are  described  as  hyperelastic. 

3. 2  Kinematics 

Let  C.l  be  the  undeformed  configuration  or  natural  state 
of  an  elastic  body.  To  be  consistent  with  the  notation  introduced  in 
Chapter  II  the  deformed  configuration  of  the  body  is  called  C.3  and 
since  the  body  is  purely  elastic  no  reference  is  made  here  to  con¬ 
figurations  C.2  and  C.4  discussed  previously  in  connection  with  elastic- 
plastic  materials. 


- 
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Let  5-j ,  anc*  £3  be  the  coordinates  of  a  material  point 
referred  to  a  curvilinear  convected  coordinate  system  [48],  that  is 
a  coordinate  system  whose  coordinate  surfaces  consist  of  the  same 
material  particles  at  all  times  during  the  deformation .  Let  the 
coordinates  of  a  particle  referred  to  a  fixed  rectangular  Cartesian 
coordinate  system  be  and  in  C.l  and  C.3  respecti vely .  The  con¬ 
vected  coordinate  system  is  defined  by 

XK  =  XR  (CrC2^3)  ,  (3.2.1) 

where  .X^  are  single  valued,  continuous  functions  of  the  £. ,  and  the 
motion  of  the  body  is  given  by 

■  X.  =  xi  (e-,  .€2>€3.t)  ,  (3.2.2) 

where  the  x^  are  continuous  and  single  valued  functions  of  £.  for 
any  fixed  t. 

If  r  is  the  position  vector  of  a  particle  P  in  C,3  referred 
to  the  origin  of  the  fixed  Cartesian  coordinate  system  and  R  is  the 
position  vector  of  the  same  particle  in  C.l,  the  displacement  vector 
is 


(3.2.3) 
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From  equations  (3.2.1)  and  (3.2.2)  it  follows  that 


R  =  R(£-j  ,52,53)  , 

r  =  r^l  *^2  ’^3*^  * 
and  u  =  u(51 ,52,53,t)  . 

If  the  vectors  e.  are  unit  vectors  along  the  coordinate 
axes  of  the  Cartesian  coordinate  system  then* 

R  =  XK  eK  , 

,  k  - 

and  .  r  -  x  e^  . 


The  local  covariant  bases  vectors  at  an  arbitrary  point  PQ  in  C.l  and 
P  in  C.3  for  this  curvilinear  coordinate  system  are  then 


and 


(*4) 

35  p, 


(^4) 

P 


(3.2.4a,b) 


*The  same__Cartesian  coordinate  system  is  used  to  describe  C.l  and  C.3 
so  that  e^  and  ex  refer  to  the  same  vectors,  however  the  lower  and 
upper  case  subscripts  are  both  used  so  as  to  be  consistent  with  the 
summation  convention. 
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respecti vely ,  where  the  comma  followed  by  the  subscript  1  denotes 
partial  differentiation  with  respect  to  £ .  and  the  subscripts  on  the 
parentheses  denote  the  points  at  which  the  enclosed  quantities  are 
evaluated. 

The  contravariant  base  vectors  in  C.l  and  C.3  are  defined 
by  the  inner  products 


0  0 

G1  .  G"  =  6nj  s 
J 

o  o 

and  g1  •  g.  s  61j  . 

Thus  the  components  of  the  metric  tensors  in  C.l  and  C.3,  referred 
to  the  convected  coordinate  system  are 


ij  =  a?1  id  G 

U  |y  1/  »  u 

3Xk  axN 

axK  axK 

J*  ac1  acJ' 

(3.2,5a,b) 

and 

ij  ar1 

9  =  k  *k  ’  gi 

axk  axk 

J  ac1  acJ” 

(3.2.6a  ,b) 

ax  axK 

so  that  the  lengths  dS  and  ds  of  an  infinitesimal  line  element  in  C.l 
and  C.3  are  given  by 


dS2  =  G  .  .dj’d^  , 

*  J 

ds2  =  g^d^d^  • 


and 


. 


35 


The  difference  of  the  squares  of  the  lengths  is 


ds2  -  dS2  =  Y.jds’dsJ  , 

where  Y^  ■  \  (9^-^)  (3.2.7) 

is  a  strain  tensor.  From  equations  (3.2.5b),  (3.2.6b),  (2.2.6a),  and 
(3.2.7)  it  follows  that 


3XK  3Xl 
3C1  3£J' 


and  from  equations  (3.2.3),  (3.2.4),  and  (3.2.7) 


t(G.+u  . ) • (G  .+u  .)-G. .] 

•  » I  j  ij  I  j 


The  material  derivative  of  this  equation  is 


(3.2.8) 


1 

2 


(v.  .+v.  •) 
j;i 


(3.2.9) 


where  v.  are  the  covariant  components  of  the  velocity  vector  in  C.3 
referred  to  the  local  basis  g1  and  a  subscript  following  a  semi-colon 
indicates  covariant  differentiation  with  respect  to  ^  in  C.3.  That  is 


5 


.  • 
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where  the  Christoffel  symbols  r?.  are  evaluated  from  the  metric  tensor 

•  J 


g 


The  components  of  y.  .  in  equation  (3.2.9)  are  the  components 

*  J 

of  the  rate  of  deformation  tensor  d.  .  referred  to  the  convected  co- 

'  J 

ordinate  system. 

Mixed  tensors  may  be  formed  from  y..  in  two  ways  by  using 

’  J 

the  metric  tensors  from  C.l  and  C„3  to  raise  one  index.  The  following 
definition  is  adopted 


•  O  I 

i  ri  k 

Y  j  - G 


-i<8VaV 


Three  independent  invariants  of  y1 .  are  of  interest.  They 

J 

are  also  the  principal  invariants  of  Green's  deformation  tensor  [49] 


3x.  3x^ 


'kl  axK  axL 


and  Finger's  deformation  tensor  [50] 


3x.  3x. 

B  .  =  — - - i 


ij  axK  axK 


and  are  the  coefficients  of  the  powers  of  n  in  the  expansion  of  the 


determinant 


fctf  nT  .xebnr  sio  92in  oJ  £  3  bne  f.D  mo- '  zio<  » •-.  Drilflfn  srtJ 


•  n  < 


- 
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|r,6’j +  +  +  Glkgkj 

■  n3  +  I^2  +  I2n  +  l3  =  0  . 


It  follows  that 


and 


where 


h  =  3  +  2A  =  e1J9ij  > 

t  _  o  1  <  k  ,  « /  k  m  k  m  \ 
I 9  =  3  +  4y  ,  +  2(y  , y  -y  Y  , ) 
2  k  '  k1  m  '  m1  k' 

1  / T  2  rimrjn  „  v 
-  y  ( 1 1  -G  G  gijgmn)  , 

I3  =  det(G1mgmj)  =  g/G  , 

G  =  detjG^)  and  g  =  det(g.j)  . 


3.3  Stress  Tensor  and  Equations  of  Motion 

The  Cauchy  stress  tensor,  which  is  symmetric,  has  components 
a.,  referred  to  a  Cartesian  coordinate  system.  Consider  the  contra- 

*  J 

variant  and  mixed  components  of  this  tensor  referred  to  the  convected 
coordinate  system  in  C.3.  From  the  transformation  law  for  second  order 
tensors  these  components  are 

Tre 

3X1  3xJ  1J 


(3.3.1) 


•  ■ 


' 


' 
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and 


rs 

It  follows  that  t  is  symmetric. 


If  T  is  the  stress  vector  acting  on  a  surface  in  C.3  with 
unit  normal  vector  n"  where 


then 


Cauchy's  first  law  applied  to  an  arbitrary  material  volume 
V  with  closed  surface  S  in  C.3  is 


T  dS  +  p3(b-a)dl/  =  0 
« 

1/ 


(3.3.2) 


where  b  is  the  body  force  vector  per  unit  mass,  a  is  the  acceleration 
vector  and  p3  the  density  in  C.3.  Equation  (3.3.2)  is  an  analogue  of 
Newton's  second  law  and  may  be  derived  from  an  energy  balance  on  the 
volume  V  using  invariance  requirements  under  rigid  body  motions  [51]. 

Applying  Greens  theorem  [52]  to  the  surface  integral  in 
equation  (3.3.2)  gives 


[xrs;r  +  P3(bs-as)]9Sdl/  =  0  • 


'  /  ■  ■ 

?j- :  .3  fi  z  o  ' •  o  >?.•  -  :  '  'J 

M  /  /  -  ”SjK  [  \J  a 

'  V, 
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Since  the  volume  V  is  arbitrary  and  the  vectors  gf5  form  a  basis  in 
C.3,  Cauchy's  first  equation  of  motion  in  terms  of  components  re¬ 
ferred  to  the  basis  g ^  in  Co 3  is 

T*s;r  +  p3^bs‘as^  =  0  ‘  (3.3.3) 

3 . 4  The  Constitutive  Equations  for  an  Hyperelastic  Material 

The  rate  at  which  the  surface  forces  do  work  on  S  plus  the 
rate  at  which  the  body  forces  do  work  on  the  mass  in  the  volume  l / 
minus  the  rate  of  increase  of  the  kinetic  energy  of  the  mass  in  l /  is 


vdS  + 


p^b«vdl/ 

d 

1/ 


^a*  vdl/ 


(3.4.1) 


s  -Tf  n  vSdS  +  f p o ( b  -a  )vsdl/. 
s  r  J  3  s  s 

5  V 

Converting  the  surface  integral  into  a  volume  integral  yields  *  after 
rearrangement 


1/ 

The  first  term  in  the  integrand  is  zero  by  equation  (3.3.3)  so  that 


H  - 


xr  vs  dl/ 
s  ;r 


TrSV  dl/  . 
s ;  r  5 


l / 


V 


" 
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and  from  the  symmetry  of  the  stress  tensor 


H  « 


rs 

T 


+ 


)  dt/  , 


or  using  equation  (3. 2.9) 


V 


Also  from  equations  (2.3.1),  (3,2.8),  and  (3.3,1)  this  may  be  written 
as 


H  = 


l / 


a.  od .  .dl/ 
ij  ij 


rs 

Materials  for  which  the  stress  t  at  a  point  is  a  single 

valued  function  of  the  deformation  at  that  point,  as  given  by  y..,  and 

*  3 

for  which 


T 


rs 


(3.4.2) 


where  E  at  a  point  is  a  continuous  single  valued  differentiable 
function  of  y-.,  are  called  hyperelastic  materials  and  the  function  E 
is  called  the  elastic  potential  or  strain  energy  function  per  unit 
mass . 


Since 


I  I 


' 


. 


•  • 
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it  follows  that 


(3.4.3) 


where  all  nine  components  of  y. .  are  taken  to  be  independent  for  the 

®  vj 

purpose  of  partial  differentiation.  From  equations  (3.4.2)  and  (3.4,3) 
the  stress  tensor  is 


rs 

T  -  p 


3E 


3  3y 


rs 


Expressed  in  terms  of  W,  the  elastic  potential  per  unit 
volume  in  C.l,  this  becomes 


(3.4.4) 


where  W  =  p-jE,  and  p-j  is  the  density  in  C.l. 

In  general  for  a  homogeneous  isotropic  hyperelastic  body, 

W  is  a  function  only  of  the  invariants  I-j  ,  I^,  I3,  and  it  can  be  shown 
that  [53] 


xrs  *  4>GrS  +  Y  brs  +  pgrs  ,  (3.4.5) 

,  _  2  3W 

~st3  81 1 


where 


* 


- 
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and 


5 


p  =  2  A 


-  3W 
3  91  ' 


krs  _  T  rrs  rmrrns 
b  =  1 1 G  -  G  G  g 


mn 


(3.4.6) 


The  constitutive  equation  (3.4.5)  referred  to  a  fixed 
Cartesian  coordinate  system  is 


where 


°ij  =  p6ij  +  <<}>+1  iT)Bi j  -  '*'Bi kBkj  > 


B 


9x.  9x. 
L 1 


ij  9Xm 


(3.4.7) 


are  the  components  of  Finger's  strain  tensor  [50]. 
For  an  incompressible  material 


and  p,  which  represents  an  hydrostatic  pressure,  cannot  be  determined 
from  equation  (3.4.6)  but  must  be  found  from  the  equilibrium  equation 
and  the  boundary  conditions  for  the  particular  problem. 

Two  particular  forms  of  the  elastic  potential  W  for  elastically 


incompressible  materials  are  that  postulated  by  Mooney  [54]  for  which 


. 


■ 
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W  =  C-j  (I-j-3)  +  C2(I2-3)  , 
and  that  for  a  neo-Hookean  material  for  which 

W  »  0(1-,  -3 )  . 

For  small  strains  the  constitutive  equation  for  a  neo-Hookean  material 
reduces  to  that  for  an  incompressible  Hookean  elastic  material  with  a 
Young's  modulus  of  6C.  A  theoretical  derivation  of  the  form  of  W  for 
a  neo-Hookean  material  is  shown  in  a  book  by  Treloar  [55]  on  the  physics 
of  rubber  elasticity. 

The  above  discussion  does  not  include  the  effects  of 
temperature  or  temperature  gradients.  A  more  general  analysis  in 
section  3.6  which  considers  thermal  effects  shows  that  for  at  least 
two  particular  types  of  loading  an  elastic  potential  W  does  exist. 

3.5  Energy  Balance  and  the  Entropy  Inequality 

The  principle  of  conservation  of  energy  states  that  in  a 
material  volume  1/  with  surface  5  in  configuration  C.3,  the  rate  of 

e  e 

change  K  of  the  kinetic  energy  plus  the  rate  of  change  U  of  the  in¬ 
ternal  energy,  must  equal  P  the  sum  of  the  rate  at  which  the  surface 
tractions  do  work  on  S  and  the  rate  at  which  the  body  forces  do  work 
on  the  mass  in  l/,  plus  the  rate  at  which  all  other  energy  such  as 
heat,  electrical,  or  chemical  energy  is  added  to  the  volume  (/, 

In  the  present  discussion  only  heat  energy  is  considered  in  the  last 


■ 


' 

tr,  %9tsm 
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term.  This  principle  may  then  be  written  as 

o  e 

K  +  U  =  P  +  d  (3.5.1) 

where  £  is  the  rate  at  which  heat  is  added  to  the  volume  l/. 

Rearranging  equation  (3.5.1)  and  recalling  the  definition 
of  H  in  equation  (3.4.1)  gives 


U  *  H  +  d  , 


or* 


Udl/ 


1/ 


a.  .d.  .dl/ 
ij  iJ 


J  qin,ds 

S 


+ 


p^h  dt/  , 


\J 


(3.5.2) 


where  U  is  the  internal  energy  per  unit  mass,  q.  are  the  components 
of  the  heat  vector  per  unit  area  at  a  point  on  S  with  unit  normal  n. , 
and  h  is  a  heat  source  per  unit  mass  in  (/. 

Conservation  of  mass  gives  that 


l / 


so  that  equation  (3.5.2)  becomes,  after  converting  the  surface  integral 


*A  superposed  dot  and  D/Dt  both  denote  material  differentiation. 


. 
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to  a  volume  integral  , 


3q . 


I / 


and  since  the  volume  V  is  arbitrary,  the  energy  equation  becomes, 
in  local  form 


(3.5.3) 


Denoting  by  S  the  entropy  per  unit  mass  and  the  temperature 
by  T,  (T>0) ,  the  following  entropy  inequality  is  postulated 


Dt 


1/ 


l / 


This  is  known  as  the  Clausius-Duhem  inequality  [56],  and  it  states  that 
the  total  rate  of  increase  of  entropy  in  the  body  is  greater  than  or 
equal  to  the  rate  at  which  entropy  flows  through  the  surface  S,  plus  the 
rate  at  which  entropy  is  produced  by  sources  within  the  body.  The  in¬ 
equality  simplifies  to 


and  using  the  energy  equation  (3.5.3)  this  becomes 


- 

' 


. 
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(3.5.4) 


Expressed  in  terms  of  A,  the  Helmholtz  free  energy  per  unit  mass  [57] 
defined  by 


(3.5.5) 


A  =  U  -  TS, 


equation  (3.5.4)  becomes 


-  p3(A  +  TS)  +  a.  .d.j  - 


(3.5.6) 


3.6  Thermoel asticity 

A  theory  of  elasticity  based  on  the  general  definition  of 
elasticity  but  which  includes  thermal  effects  is  known  as  thermo¬ 
elasticity  [58].  The  constitutive  equations  for  materials  included 
in  this  theory  may  be  written  as 


s  =  S(T>  axT5  > 


and 


. 

- 


■ 
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The  principle  of  material  indifference  gives 


and 


A  =  A(T,Ekl) 

s  =  s(t,ekl)  . 


The  material  derivative  of  the  free  energy  A  using  equation 

(2.3.1)  is 


3x.  3x 

— ! — _JL  H 
axK  9Xl  aij  • 


Substitution  into  the  entropy  inequality  (3*5*6)  gives 


"  P3^9T  +  S^.+  ^ai j”p3 


9A 


9E 


KL 


9x*  9x. 

_L 

9Xk  9Xl ' 


'i  9T 


ij 


3xi  - 


>  0  .  (3.6.1) 


Since  T  and  djj  may  be  chosen  arbitrarily,  with  the  restriction 
that  d.>  be  symmetric,  inequality  (3*6,1)  implies  that 


S  - 


9A 

9T  • 


(3*6*2) 


a 


ij 


3xi  3x.j  3A 
p3  3Xk  3X.  3Ekl 


(3.6.3) 


and 


£L 

T  9x. 


>  0  . 


48 


Considering  the  internal  energy  as  a  function  of  S  and  E 
the  entropy  inequality  (3.5.4)  implies  that 


KL 


T  = 


au 

3S  * 


and 


aij 


3xi  3xi  3U 
°3  3XK  3\  3EKL 


(3.6,4) 


Equations  (3.6.3)  and  (3.6.4)  indicate  that 


(#~>T  =  ( 

dtKL  T 


3U 


3E 


KL 


where  the  subscripts  on  the  brackets  indicate  variables  which  are 
held  constant  for  the  differentiation.  The  constitutive  equations 
(3.6.3)  and  (3.6.4)  are  valid  for  any  type  of  loading.  In  particular 
however,  equation  (3.4.4)  with  W  =  p-jA  is  a  special  case  of  equation 
(3.6.3)  for  an  isothermal  process  and  with  W  =  p^U  it  is  a  special 
case  of  equation  (3.6.4)  for  an  isentropic  process. 

3.7  Thermoelastic  Considerations  for  an  Elastic-Plastic  Material 

Consider  an  el astic-perfectly  plastic  body  with  the  con¬ 
figurations  C.l,  C.2,  and  C.3  defined  in  section  2.1.  For  convenience 
the  Piola-Ki rchoff  stress  tensor  [59]  referred  to  C.2, 


p3  3xi  3xj  ' 


(3.7.1) 


- 

. 
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and  the  stress  tensor  with  components 


S 


KL 


(3.7.2) 


( p) 

where  the  proper  orthogonal  tensor  is  defined  in  section  2.3, 
are  introduced.  The  assumption  from  classical  plasticity  that  the 
plastic  volume  change  be  zero  is  retained  so  that  in  equation  (3.7.1) 
p-|  appears  since  it  equals  p^. 

Unlike  S^,  the  components  of  remain  unchanged  during 
a  rigid  body  motion  of  C.2.  If  it  is  assumed  that  the  elastic 
moduli  and  specific  heats  are  unchanged  by  previous  plastic  defor¬ 
mation,  the  result 


SKL  ■  SKL<T’EMN}>  <3-7-3* 

is  obtained  since  for  unloading  from  the  elastic-plastic  state  the 
body  behaves  as  a  hyperelastic  material.  Furthermore  the  constitutive 
assumptions 

A  °  A(T*EKL),EKL))  (3.7.4) 

and  S  =  SCT.E^.E^b 

are  made.*  Green  and  Naghdi  [23]  and  Kestin  [60]  have  assumed  that  the 


- —  ■  -  rpy 

*A  and  S  might  depend  on  the  plastic  strain  history  as  well  as  on  . 


ksfltaido  z\ 


' 

• 

* 

* 
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entropy  is  independent  of  the  plastic  strain,  for  the  finite  and  small 
strain  elastic-plastic  theories  respectively,  however  there  seems  to  be 
no  physical  basis  for  this  assumption. 

For  elastic  unloading  of  the  body  from  C.3  the  entropy  inequality 
(3.5.6)  becomes 


-(A+TS) 


a.  .d^  q .  T 

+  IJ  _  _J_  9J. 


P3T  9x.j  — 


>  0 


(3.7,5) 


since  d* .  =0.  From  the  definitions  of  S  Q  and  d(^  it  follows  that 
ij  a$ 


J_  0  d ( e )  s  _L  5  §(e) 

P3  aij  ij  P-j  a3 


and  from  equations  (2.3.6)  and  (3.7.2)  this  becomes 


1  d(e)  .  J_s  E(e) 
P3  ij  ij  P1  KI^KL 


since  CL,  =  0  when  there  is  no  plastic  deformation  occurring.  Equation 
l\L 

(3.7.5)  may  then  be  written  as 


-(A+TS)  +  —  S  — -  >  0 


From  equation  (3.7.4)  it  follows  that  for  elastic  unloading 


A 


_  9A  *  9A  f-(e) 

'  W  T  +  TTTiT  ekl 


9E 


KL 


(3.7.6) 


so  that  the  entropy  inequality  for  unloading  is,  after  substitution 
of  (3.7.6) 


(f£+  S)T  +  (- 


3A 


3E 


1  S„,  )E,(e)  ^  3T 


TeT  +  P,  SKL,EKL  '  p-T  3X.  - 
KL  1  pi 


>  0  .  (3.7.7) 


*  -(e) 

Both  T  and  E^  may  be  chosen  arbitrarily  for  elastic  unloading  or 
reloading  before  further  yielding  so  that  equation  (3.7,7)  gives 


s  -  -  (If) 

E(e),E(P) 


(3.7.8) 


3A 


u  =  Pi(^ST>t,e(p) 


(3.7.9) 


and 


■i  3T 


PoT  8x.  - 
3  l 


>  0  . 


Equations  (3.7.8)  and  (3.7.9)  give  the  Maxwell  relation 


3S 


/  ao  \ 

"  r(P) 


9Ski 

(^) 


3Ekl  t,e 


3T  'E(e)  ’ 


and  from  equation  (3.7.3)  this  may  be  written 


,  3S  ' 

3E^p  T,E(P) 


^Kl7T  >emn^  ^  ’ 


(3.7.10) 


where  is  a  tensor  valued  function  of  the  temperature  and  the  elastic 
KL 


strain. 


. 

■  % 
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A  similar  analysis  using  the  form  (3.5.4)  of  the  entropy 
inequality  gives 


T 


■  # 


■(e) 


If  the  further  assumption  is  made  that  the  specific  heat  at  constant 
strain  [61] 


C  =  (— ) 
lE 


E(e)>E(P) 


=  t(— ) 
1 


E(e)jE(p) 


is  independent  of  previous  plastic  straining  then 


(If)  =  <f'(T,E*|e))  ,  (3.7.11) 

E(e).E(P) 

where  ^  is  a  scalar  function  of  the  temperature  and  the  elastic  strain . 
Equations  (3.7.10)  and  (3.7.11)  show  that 


S  =  S^CT.E^)  +  S(P)  (E^5 )  (3.7.12) 

fe)  (P) 

where  Sv  '  and  Sv  '  are  the  elastic  and  plastic  parts  of  the  entropy 
respectively. 

From  equations  (3.7.8)  and  (3.7.9)  the  free  energy  may  be 

written  as 

A  =  A<e>(T,E<*>)  +  A^tT.E^b 


(3.7.13) 


' 
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where  A 


is  the  elastic  free  energy  and 


A<p>  =  -  TS^  +  A^(E^) 

(P) 

is  the  plastic  free  energy  and  Av  '  is  a  function  of  the  plastic  strain 
and  possibly  the  plastic  strain  history.  From  equations  (3.7.8), 
(3.7.9),  (3.7.12),  and  (3.7.13) 


(e)  3 A ( e ) 

s  '  3T 


(3.7.14) 


c  .  3A^ 

SKL  '  P1  TTIeT 
3EKL 


(3.7.15) 


The  equation 


;(e)  .  3A  6  ;(e) 

A  ’  ZZ5T  ekl 

3EKL 


is  valid  whether  or  not  plastic  deformation  is  occurring  and  from 
equations  (3.7.14)  and  (3.7.15)  it  becomes 


1 


pi 


S  E 
5KLlKL 


-  S 


For  isothermal  deformation 


3 


- 
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*  (  6 )  *  f  6 ) 

where  Wv  '  «=  S^E^'  is  the  rate  of  increase  of  elastic  strain  energy 
per  unit  volume  in  C.l.  If  the  material  is  elastically  isotropic  S 


and  E 


(e) 


KL 

the  result 


KL 

are  coaxial  as  are  and  E^.  Thus  from  equation  (2.3.8) 


c  p  (  6  )  _  c:  p  (  6 )  -  c  p  (  ®  ) 

nLnL  a3  a3  a3  a3 


is  obtained  for  an  elastically  isotropic  material  so  that 


=  S  rE(®>  =  —  a  dj?^  . 
a3  a3  P3  ij  ij 


(3.7.16) 


The  total  rate  of  stress  work  per  unit  volume  in  C.l  is 


P1 

U  =  —  a.  .d.  .  , 
P3  1J  iJ) 


so  that  from  equation  (2.3.10) 


W  =  W 


(e) 


where 


a.  .d^r . 

ij  iJ 


(3.7.17) 


is  the  rate  of  plastic  energy  dissipation  per  unit  volume  in  C.l.  If 
the  material  is  elastically  anisotropic  then  the  rate  of  increase  of 
elastic  strain  energy  and  the  rate  of  plastic  energy  dissipation  are 


' 

• 

. 
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not  given  by  equations  (3.7.16)  and  (3.7.17). 

Elastically  isotropic  materials  only  are  considered  in  the 
remainder  of  this  thesis. 
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CHAPTER  IV 

YIELD  CONDITIONS  AND  FLOW  RULES  FOR  PLASTIC 
DEFORMATION  WITH  FINITE  ELASTIC  STRAINS 

4.1  Introduction 

A  flow  rule  is  developed  in  section  4.2  which  is  based  on  the 
assumption  that  the  materials  to  be  considered  satisfy  Drucker's 
postulate  [30].  The  further  assumptions  of  the  existence  of  a  yield 
function,  no  work-hardening,  incompressibility,*  elastic  isotropy,  and 
isothermal  deformation  are  made.  In  section  4.3  the  theory  is  limited 
to  the  discussion  of  materials  whose  elastic  behavior  is  neo-Hookean 
and  two  yield  conditions  which  are  related  to  the  classical  von  Mises 
yield  condition  are  discussed.  In  section  (4.4)  a  simplification  which 
results  from  the  further  assumption  of  plastic  isotropy  is  indicated. 

4.2  A  Plastic  Flow  Rule 

In  the  classical  elastic-plastic  theory  for  perfectly  plastic 
solids  it  is  assumed  that  there  exists  a  yield  condition 

f(a.j,T)  =  0  (4.2.1) 

with  the  elastic  domain  defined  by  f(o^.,T)  <  0  . 


*The  theory  may  be  generalized  to  include  finite  elastic  volume  change 
[76]. 


\:z  :  v  ±_ 
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If  a  material  is  plastically  isotropic  the  form  of  the  function  f  at 
a  point  must  be  independent  of  any  rigid  body  rotation  of  the  body 
and  f  is  then  a  function  of  the  invariants  of  the  stress  tensor. 

However  for  a  plastically  anisotropic  material  there  are 
two  situations  for  which  the  form  of  the  yield  function  f(a.  .  ,T)  is 

*  J 

not  satisfactory.  If  the  elastic  strains  just  prior  to  yield  are 
small  as  in  classical  elasticity  but  the  rotations  are  large,  as  in 
the  problem  of  the  elastica  [6],  or  if  the  material  undergoes  finite 
elastic  strains  and  rotations,  the  form  of  the  yield  function  at  a 
point  in  a  plastically  anisotropic  material  must  depend  on  the  rotation 
of  the  material  element  at  that  point.  The  yield  condition  is  then  of 
the  form 


f(aKL,T)  =  0  (4.2.2) 

proposed  by  Green  and  Naghdi  [23],  The  stress  cy,  is  given  by 

aKL  =  RiKRjLai j  ’ 

where  R.  ^  is  a  proper  orthogonal  tensor  representing  a  pure  rotation 
and  is  obtained  from  the  polar  decomposition  of  the  total  deformation 
gradient.  That  is 


FiL  RiKUKL  * 


iwsna  *d  bsaoqoiq 
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where  is  a  symmetric  positive  definite  tensor  representing  a  pure 
stretch  deformation.  Thus  ^re  the  components  of  the  stress  tensor 
at  a  point  referred  to  a  Cartesian  coordinate  system  which  has  under¬ 
gone  the  same  rotation  as  the  material  element  at  that  point. 

If  large  plastic  deformations  occur  the  yield  condition  for 
a  plastically  anisotropic  material  is  of  the  form  (4.2.2)  regardless 
of  the  magnitude  of  the  elastic  strains. 

(e) 

The  deformation  rate  tensors  d..,  d.  and  d* .  may  also  be 

*  J  *  J  *  J 

referred  to  coordinates  which  have  undergone  the  same  rotations  as 
the  material  elements  so  that 


[dKL’dKL^  ,dKL^ 


RiKRjlTdi j  ,di j  ^  ’dij^ 


and  since  is  an  orthogonal  tensor  it  follows  that 


°KL^dKL,dKL^,di<L^  " 


a.  .[d.  .,d!^,d*.] 
i J L  i J  *  lj  i J J 


Consider  a  body  which  is  both  elastically  and  plastically 
incompressible  and  elastically  isotropic.  Assume  that  a  yield  function 
f(al^,T)  exists  and  suppose  that  at  time  t°  the  body  is  in  equilibrium 
at  the  temperature  T  and  homogeneous  stress  which  lies  inside  the 
yield  surface  but  sufficiently  close  to  the  yield  surface  that  the  ad¬ 
dition  of  the  stress  necessary  to  move  the  stress  point  to  the  yield 
surface  results  only  in  infinitesimal  additional  elastic  strains.  Let 


,.5.t  fciHK .iWJav'  6fl  t.uca  ns  »  l~)  6  oJ  Hi*'  «  6  *6 


‘  ",  •  [t|b. ' ,;  f4«  i  fb  r  r°  *  j$  '  J>  ’•  JX"  :  JX! 

. 

. 


59 


an  external  agency  apply  isothermally  an  additional  stress  to  the  body 
so  as  to  move  the  stress  point,  as  shown  in  Figure  4.1,  to  a ^  at  time  t*, 
and  then  along  the  yield  surface  to  the  point  a |V  at  time  t' 1  =  t*  +  6t, 
and  finally  back  to  the  original  stress  a°L  at  time  t.  In  the  special 
loading  case  considered  here  the  restriction  is  made 


Fi gure  4. 1 

A  Diagrammatic  Representation  of  the 
Isotropic  Yield  Surface  in  Stress  Space 

that  the  plastic  strains  occurring  during  the  time  interval  [t*,t''] 
are  infinitesimal  and  since  the  elastic  strains  are  similarly  re¬ 
stricted,  changes  in  the  configuration  of  the  body  during  the  cycle 
may  be  neglected. 

According  to  Drucker's  postulate  the  work  done  by  the  external 
agency  during  any  complete  cycle  of  loading  and  unloading  is  non-negative. 
The  total  work  done  per  unit  volume  by  both  the  original  stress  and 
the  external  agency  is  given  by 


. 

- 

•  <  ;  d  •• 
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W 


aKLdKLdT 


.0 


t 

j  aKL^dKL^+dKL^  dTo 

u.0 


Since  the  plastic  deformation  occurs  only  from  time  t*  to 
t*+6t,  and  d*  =  0  at  all  other  times,  the  expression  for  the  total 
work  done  is 

t*+St 

W  s  ^  aKldKL  d  +  j  aKLdKLdT* 

t* 

where  0  denotes  integration  around  the  complete  cycle  of  loading 

T 

and  unloading. 

Since  the  material  is  elastically  isotropic  and  incompressible 
it  follows  from  the  discussion  in  section  3,7  that 


d(e) 
KL  KL 


dx 


and  thus 


t*+6t 


°KLdKLdT' 


Furthermore  since  the  changes  in  configuration  of  the  body  may  be  neg¬ 
lected  it  follows  that  the  work  done  by  the  initial  stress  is  given 


. 

'  "  si 


- 
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by 


t*+6t 

- 

e 


0  H 

GKLaKL 


di  . 


Thus  the  net  work  done  by  the  external  agency  during  the  complete 
cycle  is 

t*+fit 

W  -  W°  ^  j  ■ aKL”aKL^dKLdx  * 

t* 


From  Drucker's  postulate 


so  that 


consequently 


W  -  W°  >  0 

1 i m  W-W°  n 
6t^0  fit  - 


^°KL“aKL^  dKL  —  °  s  (4,2,3) 


where  is  a  stress  point  on  the  yield  surface  associated  with  the 
plastic  rate  of  deformation  tensor  dj*^  and  is  any  stress  point 
near  aj^»  on  or  just  inside  the  yield  surface.  Equation  (4,2,3)  implies 
that  the  vector  representing  d*^  in  nine  dimensional  stress  space  is 
normal  to  the  yield  surface  so  that 


' 


.  s"l? 
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=  A 


9f 


9a 


KL 


(4.2.4) 


where  A  is  a  non-negative  factor  of  proportionality.  For  a  pi asti cal ly 
isotropic  material  this  flow  rule  may  be  written  as 

di  j  =  X  3o77  •  <4'2'5> 

Equations  (4.2.4)  and  (4.2.5)  and  the  assumption  of  plastic 
incompressibility  imply  that  the  yield  function  must  be  independent 
of  the  hydrostatic  pressure. 

4.3  Two  Yield  Conditions  and  Associated  Flow  Rules 

The  remainder  of  this  thesis  is  restricted  to  the  discussion 
of  elastic  perfectly-pl as  tic  materials  whose  elastic  behavior  is  neo- 
Hookean.  That  is,  there  exits  an  elastic  potential  given  by 


W  *  c(l{e)-  3)  ,  (4.3.1) 

where  the  invariant  lje*,  analogous  to  I,  defined  in  section  3.2,  is 
the  first  principal  invariant  of  the  tensor 


r(e) 

a6 


so  that 


' 


■ 


.  t  I , 
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W  =  2CE*e-  . 

aa 

The  von  Mises  yield  condition  is 

7  aijaij  -  k2  ■  o  >  <4-3-2> 

where  k  is  the  yield  stress  in  pure  shear.  In  the  classical  small 
strain  elastic-plastic  theory  the  condition  (4.3,2)  is  equivalent  to 
the  condition  that  yielding  occurs  when  the  elastic  shear  strain 
energy  reaches  a  critical  value.  The  maximum  shear  strain  energy 
yield  condition  is  not  equivalent  to  (4.3,2)  for  a  material  which  may 
undergo  finite  elastic  strains  and  the  two  yield  conditions  are  hence¬ 
forth  referred  to  as  the  maximum  shear  strain  energy  yield  condition 
and  the  von  Mises  yield  condition  respecti vely . 

The  flow  rule  for  a  plastically  isotropic  material  which 
obeys  the  von  Mises  yield  condition  is,  from  equations  (4,2,5)  and 
(4.3.2), 


dtj  -  Aal j  .  (4.3.3) 

For  the  maximum  shear  strain  energy  yield  condition  yielding 
occurs  when  the  elastic  strain  energy  (equal  to  the  shear  strain  energy 
since  elastically  incompressible  materials  only  are  considered)  reaches 


a  critical  value  W*. 


. 


■ 
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From  equation  (3.4.5)  for  a  neo-Hookean  material 


a!.  =  2C[F^eM.e)  -  i  I,(e)6..]  . 

ij  la  ja  31  ijJ 


(4.3.4) 


Using  this  constitutive  equation  and  equation  (4.3.1),  the  elastic 
shear  strain  energy  may  be  expressed  as  a  function  of  the  invariants 
of  the  stress  deviator  tensor.  Define  two  new  stress  invariants* 


K9  -  -Kf  o !  .  a  i  . 

2  8C2  ij  ij 


1 


4C 


2 


and 


K 


1 


i  _.i  i 

a.  .a.,  a,  . 


1  ,, 


3  24q3  ij  jk  ki  8C3  3 


j:  . 


It  may  be  shown  that 


K, 


1  T(e)' 
3  A1 


(e) 


(4.3.5) 


and 


K  _  i  1  T(e)T(e)  _2_  T  (e)’ 
r3  1  “  3  A1  a2  27  1 


(4.3.6) 


Using  equation  (4.3,5)  to  eliminate  from  equation  (4.3.6) 
and  rearranging  gives 


3 

K2I{e)  -  7  I|e)  +  3(1 -K3)  =  0 


*K2  and  are  called  stress  invariants  even  though  they  are  dimension¬ 
less.  since  they  are  derived  from  the  stress  deviator  invariants  J9 
and  ji  .  2 


■* 
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or  K2(£  +3)  -  5  (£  +  3)3  +  3(1-K3)  =  0  . 

Consequently  the  maximum  shear  strain  energy  yield  condition  for  a 
neo-Hookean  material  is 

f  =  K2(^+  3)  -  1  (^-  +  3)3  +  3(1-K3)  =  0  (4.3.7) 

and  it  may  be  verified  that  f  <  0  when  W  <  W*  . 

Also  if  yielding  occurs  after  infinitesimal  elastic  strains 
then  from  equation  (4,3,4) 


ai  j  =  4C  eij  ’ 


W  *  4r  a!  .-a' 


8C  i j  ij  * 


and  W*  «  C  .  Therefore  the  yield  condition  (4.3,7)  becomes 

3K2  -  3  ^  +  0(~)  =  0 

c 

2 

so  that  if  terms  0(^p— )  are  neglected  this  yield  condition  reduces  to 

C 

the  von  Mises  yield  condition. 

The  yield  function  f  for  the  maximum  shear  strain  energy 
yield  condition  is  a  function  of  the  two  stress  invariants  K2  and  K^, 
and  since  the  material  is  incompressible,  the  associated  flow  rule 


' 
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(4.2.5)  for  this  yield  condition  is 


dij  A[3K 


3f  ^ij. 


3f 


2  4C‘ 


3K 


(■ 


1 


3  8C 


3  ai kakj 


K2 

3C 


-  X[(4^+  3) 


1J 


4C' 


-  ( 


3ai  kak,i 

8C3 


K„ 

rV]  • 


(4.3.8) 


4.4  A  Further  Simplification  Due  to  Isotropy 

Using  equation  (2.3.11)  the  flow  rule  for  a  plastically  iso¬ 
tropic  material  is 


1  [F(e 

2  ia 


a3  Sj 


F(e)d(P)F(e 

j6  a  6  ai 


-1 


]  -  A 


3f 


9aij 


A  significant  simplification  of  this  flow  rule  results  if  the  elastic 
deformation  gradient  tensor  is  expressed  in  the  polar  form 


p  (  6  )  -  W  (  ®  )  n  (  ®  ) 

Fia  ~  Vij  ja  * 


(e) 

where  R.  J  is  the  previously  defined  proper  orthogonal  tensor  repre- 
Ja 

(e) 

senting  the  rigid  body  rotation  of  Co 3  relative  to  C.2  and  .  is  a 
symmetric  positive  definite  tensor  such  that 


V  ( e )  v  ( e ) 

vik  vkj 


F  ^ e )  F  ^ e ) 
la  ja 


(4.4.1) 


With  the  use  of  equation  (4.4.1)  the  constitutive  equation 
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(4.3.4)  becomes 


o\.  =  2C -  i  (V(e)v(e))6>  j 
ij  L  lk  kj  3  v  mk  mk  '  ijJ 


(4.4.2) 


and  the  expression  for  d*.  becomes 

*  J 


dtj  -  1  CVik)dkm)v>^)"1  +  Vjk)dkMt)’1]  > 


(4.4.3) 


where 


d<P>  =  R^R^d^ 
km  ka  m3  a3 


(4.4.4) 


■  (e) 


From  the  constitutive  equation  (4.4.2)  the  tensors  a!,  and 

'  J 


Vi.  are  coaxial  and  since  the  materials  considered  are  plastically 
isotropic,  f  is  a  function  of  the  stress  invariants  and  the  tensors 
a!,  and  d*.  are  also  coaxial.  It  may  thus  be  deduced  from  equation 

'  J  *  vJ 

( p  \ 

(4.4.3)  (see  Appendix  A)  that  the  tensors  d* .  and  d'.'  are  also  co- 

1  J  *  J 

axial  so  that  equation  (4.4.3)  reduces  to 


d*. 

ij 


,<p> 

ij 


(4.4.5) 


for  a  plastically  isotropic  material. 
Furthermore  since 


F.,  -  V(eMeMf^ 

lK  ij  ja  aK 


,(P)  -  1  /r(P)F(P)‘’  r(P)r(P)'\ 

a6  2  1  ak  rKB  BK  Ka  '  ’ 


and 
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it  follows  from  equation  (4.4.3)  and  (4.4.4)  that 


—  (-V(e)  V(e) 
2  v i k  vkj 


.  V(e)"Ve) 

Vjk  vki 


v(e) 

vik 


-1 


F  F"^V^e^ 


+  V 


(e) 

jk 


-1 


F  F~V^) 
rkLrLmmi  ‘ 


(4.4.6) 
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CHAPTER  V 

SIMPLE  SHEAR  PROBLEM 

5.1  Introduction 

Simple  shear  of  a  cuboid  is  considered  in  this  chapter. 
The  sides  in  C.l  are  parallel  to  the  axes  of  a  fixed  Cartesian  co¬ 
ordinate  system.  After  deformation  the  position  of  a  particle 
originally  at  is  given  by 


X1  =  xi  +  kx2  , 


x 


2 


3 


and  x^  =  X^  . 

This  is  a  homogeneous  isochoric  deformation  and  the  parameter 
K  =  tan  0,  where  6  is  the  angle  of  shear  as  shown  in  Figure  5.1,  is 
used  as  a  time  scale.  The  deformation  is  assumed  to  occur  isothermal ly. 

The  total  deformation  gradient  and  its  total  derivative  with 
respect  to  the  time  like  parameter  K  are 


"l 

K 

o' 

0 

1 

0" 

1 — 1 
~n 

1 _ 1 

11 

0 

1 

0 

« 

and  [F-k]  = 

0 

0 

0 

l 

_0 

0 

1_ 

0 

0 

0 
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The  deformation  is  completely  specified  by  the  parameter  K. 

/ 

The  material  is  assumed  to  be  el asti c-perfectly  plastic 

and  to  exhibit  neo-Hookean  behavior  for  loading  up  to  the  yield  point 

and  for  unloading  from  an  elastic-plastic  state.  Both  the  von  Mises 

and  the  maximum  shear  strain  energy  yield  conditions  are  considered. 

The  elastic-plastic  simple  shear  problem  consists  of  finding 

(e) 

the  stress  deviator  tensor  and  the  elastic  stretch  tensor  V..  as 

*  0 

functions  of  the  parameter  K,  Complete  solutions  are  found  numerically 
for  both  yield  conditions. 
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5.2  Solution  Prior  to  Yielding 

Until  the  yield  condition  is  reached  the  cuboid  remains 
elastic  so  that  upon  removal  of  the  stress  it  returns  to  configuration 
C.l.  The  elastic  solution  was  obtained  by  Rivlin  [63]  for  both  com¬ 
pressible  and  incompressible  hyperelastic  materials. 

From  equations  (4.4.1),  (4.4.2),  and  (5.1.1a) 


so  that 


1+K2 

K 

0 

"l 

0 

o" 

C°i j]  =  2C 

K 

1 

0 

K2 

-  2C(1  + 
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_0 
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_0 

0 

1_ 

’ll 


4C  yZ 
3  K 


_  2C  „2 
22  °33  ~  "  3 


o\z  =  2CK  , 


(5.2. 1  a  ,b  ,c) 


and  all  other  stress  deviator  components  are  zero. 

If  the  stress  on  one  of  the  faces  of  the  cuboid  is  specified 
then  the  components  of  the  stress  tensor  may  be  determined  from  equations 
(5 . 2 . 1  a ,b  ,c) . 

Consider  first  the  von  Mises  yield  condition  (4.3.2)  which 
using  equations  (5.2.1a,b,c)  becomes 
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K4  +  3K2  -  |  (jr)  =  0  . 

Solving  for  the  value  of  K  at  yield  gives 

3  1  /  77?  V2 

Ky  =  (“  f  +  I  +  3  (£)  )  •  (5.2.2) 

For  the  maximum  shear  strain  energy  yield  condition  K 
must  be  determined  as  a  function  of  .  From  equation  (5.1.1a)  it 
follows  that 


I 


(e) 

1 


2 

=  3  +  r 


5 


consequently  at  yield 


or 


(5.2.3) 


A  basis  for  comparison  is  needed  for  the  solutions  with 
the  two  yield  conditions  and  the  yield  stress  in  simple  tension  is 
chosen  arbitrarily  to  be  this  basis. 

For  this  purpose  consider  the  uniform  extension  [64]  of  a 
block  of  an  incompressible  neo-Hookean  material  due  to  a  force  parallel 
to  the  x-j  direction.  The  deformation  is  specified  by 
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and 


x-j  =  AX-| 

x  .  _Lv 

2  /T  2  ' 

X  =  J.V  ' 

3  /T3’ 


so  that  the  deformation  gradient  has  components 


A  0  0 


]//R  0 
0  1//T 


It  follows  from  equation  (4.3.4)  that 


[aij] 


2CA2+p  0  0 

0  f-  +p  0 

o  0  T+P 


(5.2.4) 


If  the  faces  parallel  to  the  direction  of  extension  are  stress  free 
then 


and 


an  =  2C(A2  -  1)  . 


(5.2.5) 


o‘’*i  tb 
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The  yield  stress  Y  in  simple  tension  for  a  von  Mises  material 
is  given  by 


y=  -7k  .  (5.2.6) 

Consequently  if  the  material  is  neo-Hookean  up  to  yielding,  it  will 

yield  in  simple  tension  when  the  extension  is  the  solution  A  of  the 

Y 

equation 

Ay  ‘  T  (ir}  Ay  -  1  =  0  >  (5.2.7) 

which  results  from  equations  (5.2.5)  and  (5.2.6). 

Using  equation  (5.2.4)  the  invariant  lje^is  found  to  be 

lje)=A2+f  (5.2.8) 

for  uniform  extension.  Thus  if  a  material  whose  yielding  is  governed 

by  the  maximum  shear  strain  energy  yield  condition  is  to  yield  in 

simple  tension  at  the  same  stress  and  therefore  at  the  extension  A 

Y 

W* 

found  from  equation  (5.2.7),  the  material  constant  (--)  must  satisfy 
the  equation 


(5.2.9) 


' 
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In  the  remainder  of  this  chapter  and  in  Chapter  VI,  all 
the  stresses  are  non-dimensional i zed  by  division  by  the  shear  yield 
stress  k.  The  values  of  k  used  for  the  materials  obeying  the  maximum 
strain  energy  yield  condition  are  those  determined  from  equations 
(5.2.7)  and  (5.2.9)  for  von  Mises  materials  which  have  the  same  yield 
stress  in  simple  tension. 


5.3  Material  Obeying  the  von  Mises  Yield  Condition 

The  elastic-perfectly  plastic  solution  is  first  obtained 
for  a  material  which  is  neo-Hookean  up  to  the  yield  point  and  during 
unloading  from  the  yield  point,  which  is  determined  by  the  von  Mises 
yield  condition. 

Using  equations  (4.4.2)  and  (4.4.6)  the  flow  rule  (4.3.3) 

becomes 


y(e)  +  (V^  V^)  =  FF~^e) 


T 

) 


4A 

"  TkTcJ 


(5.3.1) 


which  is  written  in  matrix  notation  with  the  superscript  T  denoting 
the  transpose,  I  the  identity  matrix,  and  tr(.)  the  trace  of  the 
argument  matrix.  The  non-negative  parameter  A  in  equation  (5.3.1) 
differs  from  that  in  equation  (4.3.3)  by  a  factor  k. 

The  von  Mises  yield  condition  after  substitution  of 


equation  (4.4.2)  becomes 
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tr[ - - 

(k/C) 


2 


+ 


)I}] 


2  =  0 


and  taking  the  total  derivative  with  respect  to  the  time  like  para¬ 
meter  K  gives 


w  Ctr(v(e)4)  -  ltr2(v(e)2)]  =  0 

or  tr(V(e^3V^)  -  A  tr(V^e'2)tr(V*e^V^eh  =  0  .  (5.3.2) 

•V  t  *V  'V 

The  elastic  solution  (5.2.1)  is  valid  until  K  reaches  K 

as  determined  from  equation  (5.2.2).  For  K  =  K,  the  matrix 

y  ~ 

and  the  components  of  the  stress  deviator  are  known.  For  K  > 

(e) 

the  components  of  V  are  determined  from  the  system  of  non-linear 

ordinary  differential  equations  (5.3.1)  and  (5.3.2)  using  the  value 
(e) 

of  Vv  '  at  K  =  K  as  the  initial  condition.  The  numerical  method 
used  to  obtain  a  solution  is  described  in  section  (5.5). 

5.4  Material  Obeying  the  Maximum  Shear  Strain  Energy  Yield  Condition 

The  solution  for  the  simple  shear  of  a  cuboid  of  a  neo- 
Hookean  el asti c-perfectly  plastic  material  which  obeys  the  maximum 
shear  strain  energy  yield  condition  is  obtained  in  a  manner  similar 
to  that  discussed  in  the  previous  section  for  the  von  Mises  yield 
condition. 

Substitution  of  equation  (4.4.2)  into  the  expressions  (4.3.5) 


■ 


77 


and  (4.3.6)  for  the  stress  invariants  and  gives 


K2  =  1  [tr(v(e)4)  -  1  tr2(v(e)2)] 


(5.4.1) 


and 


6  4  2 

K,  =  1  [tr(V(e)  )  -  tr (V(e)  )tr(V(e)  ) 


+  I  tr2  (V(e)  )]  . 


(5.4.2) 


Using  equations  (4.4.2),  (4.4.5),  and  (4.4.6)  the  flow  rule  (4.3.8) 
becomes 


V  ( e )  +  (V^  V^e  ^ )  =  FF"Me^  +  (V^  FF~Meb 


/  <v 


/  *\/  *V 


-  x(fe{(^  +  3)[V(e)2  -  \  tr(V(e)2)I]  -  3[V(e)2  -  \  tr(V(e)2)l]2 

^  ^  ~  J  ^  ^  ~  ^  ^ 


2  2 

+  tr([V(e)  -  1  tr(V(e)  )I]2)} 


(5.4.3) 


and  the  maximum  shear  strain  energy  yield  condition  (4.3.7)  is,  after 
substitution  of  equations  (5.4.1)  and  (5.4.2) 


1  (!£-  +  3)[tr(V(e)4)  -  ]•  tr2(V(e)2)]  -  I  (£  +  3)3  +  3 


[tr(V(e)6)  -  tr(V(e)4)tr(V(e)2)  +  f  tr3(V(e)2)]  =  0  . 


' 

■ 
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Taking  the  total  derivative  with  respect  to  K  gives 

(!£■+  3)[tr(v(e)3y(e))  - 1  tr(v(e)2)tr(v(e)y(e))] 

+  [tr(v(e)4)  -  I  tr2(v(e)2)]tr(y(e)v(e)) 

5  3  2 

-  3tr(v^e*  y(eb  +  2tr(v^e^  y^)tr(y^  )  =  o  .  (5.4,4) 

Equations  (5.4.3)  and  (5.4.4)  are  solved  numerically  for 

(  0  \  ( S  ) 

Vv  '  as  a  function  of  K  using  the  value  of  Vv  '  at  yield  as  the  initial 
condition. 

5.5  Numerical  Solution 

The  two  systems  of  non-linear  ordinary  differential  equations 
(5.3.1),  (5.3.2)  and  (5.4.3),  (5.4.4)  are  sufficiently  complex  that 
closed  form  solutions  do  not  appear  possible.  Instead  a  numerical 
technique,  which  is  described  in  this  section,  was  used  to  obtain 
the  solutions. 

The  following  discussion  is  for  the  el asti c-perfectly  plastic 

solution  for  the  von  Mises  material  and  the  same  method  was  used  for 

the  material  with  the  maximum  shear  strain  energy  yield  condition. 

Consider  the  matrix  equation  (5.3.1).  Since  both  sides  are 

symmetric  3x3  matrices  it  represents  six  independent  equations  with 

(e) 

seven  unknowns;  the  six  independent  components  of  Vv  '  and  the  para¬ 
meter  X.  The  yield  condition  provides  the  additional  equation  (5.3.2) 


Jedl  xsi'q  »d  ^1  ii  2  i  « (€.^.2,  bn6  .  .  «  •  -• ;‘ ■ 

. 

' 
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for  the  determination  of  A. 

The  Runge-Kutta  method  [65],  generalized  for  a  system  of 
simultaneous  first  order  differential  equations  was  used  to  solve 
equations  (5.3.1)  and  (5.3.2).  Consider  a  system  of  m  differential 
equati ons . 


£y(1)  =  g(1)  (t.y<V2) . y(m))  , 


^y(m)  -  g(m)(t,y(1>,y<2> . ,w)  . 


(5.5.1) 


For  an  integration  step  size  At,  the  standard  fourth  order  Runge-Kutta 
formula  at  the  nth  integration  step  is 


y<]j  -  y'i}  +  f  (h|i}  +  2h<1>  +  2h<1>  +  hj1'))  +  o(At5) 


h(D  =  q(D  (t  yO)  y(2)  Jm)) 

nl  9  [Zn’yn  ,yn  •  *  *  *  ,yn  ;  » 


where 


- 


■ 
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h(i)  _  a(i)/t  +  At  (1)  At  .(1)  (2)  At  .  (2)  (m)  At  h (m) x 

h2  'g  (tn+T’yn  +  Thl  ’yn  T  hl  •  *  *  *  ’yn  "2~  hl  )s 

h(i)_  a(1)/t  .  At  (1)+  At  .(1)  (2)  At  .(2)  (m)  At  .(mh 

"3  "g  (tn+T’yn  +Th2  >yn  T  h2  •  *  *  *  >yn  T  h2  }’ 


and 


hj^  =  (tn+At,y|fj1  ^  +  Ath^,  y^  +  Ath^  ,  .  ..,y^  +  Ath^) 


The  derivation  [65]  of  these  formulae  is  lengthy  and  is  not  included 
here. 

The  system  of  differential  equations  (5.3.1)  and  (5.3.2) 

•(e) 

cannot  be  put  in  the  form  (5.5.1)  since  Vv  y  cannot  be  written  ex- 

(e) 

plicitly  as  a  function  of  K  and  V  .  For  given  values  of  K  and 

( e ) 

Vv  '  however,  it  is  possible  to  determine  the  value  of  the  matrix 
from  equations  (5.3.1)  and  (5.3.2)  by  guessing  a  value  for  A 
and  solving  equation  (5.3.1)  as  a  system  of  algebraic  equations  for 

•(e)  '(e) 

the  components  of  Vv  .  The  resulting  values  of  Vv  '  and  the  given 

values  of  are  then  substituted  into  equation  (5.3.2).  If  the 

latter  equation  is  not  satisfied,  then  the  value  of  A  is  varied  until 

•  (e) 

the  resulting  matrix  Vv  '  does  satisfy  equation  (5.3.2).  It  is  thus 

*  ( G ) 

possible  to  determine  the  value  of  the  matrix  '  as  a  function  of 
V(e)  and  K  so  that  from  a  numerical  point  of  view  the  system  of 
equations  (5.3.1)  and  (5.3.2)  is  equivalent  to  a  system  of  the  form 


ill  9nim-i9jsb  Ol  si  dr 


. 
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V(e)  =  G(K,V(e)) 

even  though  the  explicit  form  of  the  matrix  function  G  cannot  be  found. 

Thus  the  fourth  order  Runge-Kutta  formula  for  this  problem, 
in  matrix  notation  at  the  ntn  integration  step  is 

+f(Bl  +  2”2  +  2«3  +  S4>’  (5-5‘2) 

where  H,  =  G(K .V^) 

H2  =  G(K  +  y  ,  V*e)+  y  H1 )  . 

y3  =  S(K  +  f  .  v(ne)  +  f  Hz)  . 

H4  =  G(K+AK ,  v<e)+  AK  H3)  , 

and  AK  is  the  integration  stepsize. 

A  computer  program  was  written  employing  the  above  method 

to  solve  numerically  the  systems  of  equations  for  both  yield  conditions. 

Solutions  were  found  for  various  values  of  k/C  and  the  corresponding 
u* 

values  of  ,  with  total  integration  intervals  extending  to  values  of 
K  as  high  as  two. 

The  results  are  shown  graphically  in  Figures  5.3  to  5.8. 
Estimates  of  the  truncation  error  due  to  neglect  of  terms 


(  H  »A  *  (V  ,*♦*$*  *  fH 


" 

. 
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c 

0 ( AK  )  in  equation  (5.5.2)  are  available  [66],  but  since  they  are  for 
the  most  general  case,  that  is  a  system  such  as  (5.5.1),  they  tend  to 
be  pessimistic  and  rather  difficult  to  apply  to  a  complex  system  such 
as  (5.3.1)  and  (5.3.2).  An  indication  of  the  truncation  error  may 
be  obtained  however  by  changing  the  integration  step  size  AK. 

All  the  simple  shear  solutions  found  were  obtained  using 
step  sizes  of  0.01  and  0.005.  A  comparison  of  the  corresponding  re¬ 
sults  shows  complete  agreement  to  four  significant  figures  in  the 
(e) 

values  of  Vv  '  and  the  maximum  difference  observed  in  the  values  of 

-4 

the  non-dimensional i zed  stress  deviator  is  10  . 

In  an  attempt  to  reduce  the  possibility  of  round  off  error 
affecting  the  results,  double  precision  (16  significant  figures)  was 
used  for  all  calculations  and  it  is  assumed  that  the  effects  of  round 
off  error  are  negligible. 


5.6  Discussion 

The  results  indicate  that  as  k/C  becomes  smaller  the  normal 
stresses  become  smaller  and  are  negligible  for  k/C  «  1.  Also  for 
k/C  «  1  the  results  approach  those  obtained  from  the  classical  elastic- 
perfectly  plastic  theory. 

For  the  von  Mises  yield  condition  the  stress  deviator  com¬ 
ponent  approaches  zero  for  values  of  K  large  compared  with 
and  if  a ^  =  0  then  =  -  a ^  =  °\]  for  large  K.  The  maximum 
shear  strain  energy  yield  condition  gives  quite  different  normal 
stress  effects  during  elastic-plastic  flow.  The  relation  =  033 


. 
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which  holds  for  elastic  deformation  prior  to  yielding  also  holds  for 
elastic-plastic  deformation  for  all  K  and  if  a ^  =  0  then  a ^  =  0 
and  a-|  i  =  -  3a^  • 

One  further  result  is  of  interest.  Figure  5.2  which  for 
simplicity  considers  only  the  x-j  and  x^  directions,  shows  a  typical 
stress  deviator  system  in  C.3  during  elastic-plastic  flow  for  either 
the  von  Mises  or  maximum  shear  strain  energy  yield  conditions.  The 
figure  shows  the  configuration  C.2  which  results  from  irrotational 


Figure  5.2 

Direction  of  Maximum  Shear  Stress  and  the  Orientation 
of  C.2  Resulting  from  Irrotational  Elastic  Unloading 


elastic  unloading  from  C.3  with  a  being  the  angle  through  which  lines 
of  constant  rotate  during  unloading  from  C.3  to  C.2.  The  angle  <j> 
indicates  the  direction  of  maximum  shearing  stress  in  C.3.  From  the 


' 


as^rf  br  r  -  t  nh.1  r  u  , )  rrt  e  oWafifs 


graphs  of  i  al .  versus  K  it  is  seen  that  <j>  approaches  a  constant  val 

K  1 J 

for  large  K  since  the  stress  deviator  components  approach  constant 
values.  It  is  found  however  that  a  approaches  4>  as  K  becomes  large 
compared  to  K  so  that  C.2  tends  to  align  itself  so  that  the  axes 

%7 

ox*  are  in  the  directions  of  the  maximum  shearing  stress. 


; 
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Figure  5.3  von  Mises  Yield  Condition,  o\  ./k  versus 
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Figure  5.4  von  Mises  Yield  Condition,  a  I  - /k  versus 
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Figure  5.5  von  Mises  Yield  Condition,  a',  ./k  versus 
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Figure  5.6  Maximum  Shear  Strain  Energy  Yield  Condition, 
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Figure  5.7  Maximum  Shear  Strain  Energy  Yield  Condition, 
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Figure  5.8  Maximum  Shear  Strain  Energy  Yield  Condition, 
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CHAPTER  VI 


TORSION  OF  A  CYLINDER 


6.1  Introduction 


The  problem  of  the  torsion  of  a  cylinder  is  closely  related 


kinematically  with  the  simple  shear  problem  of  Chapter  V.  There  is 
the  added  complication  however  that  the  deformation  is  non-homogeneous 
so  that  use  of  the  equilibrium  equation  is  necessary.  As  for  the 
simple  shear  problem  the  deformation,  which  is  assumed  to  occur  iso¬ 
thermal  ly,  is  specified  and  the  stresses  required  to  maintain  the 
deformation  are  found  using  the  constitutive  equations  developed 
previously  for  finite  strain  el asti c-perfectly  plastic  problems. 


Consider  a  cylinder  of  length  l  and  radius  a  as  shown  in 


Figure  6.1 . 


2a 


Figure  6.1 


Torsion  of  a  Cylinder 


.V  igtfqbrtO  tftef 


92 


The  axis  of  the  fixed  Cartesian  coordinate  system  coin¬ 


cides  with  the  axis  of  the  cylinder  and  the  coordinate  system  (rs0,z) 
is  a  cylindrical  polar  coordinate  system  as  shown. 


The  deformation  is  such  that  planes  perpendi cul ar  to  the 


axis  of  revolution  remain  plane  and  all  material  particles  in  any 
such  plane  move  in  concentric  circles  with  centre  on  the  axis  of 
revolution.  This  is  an  isochoric  non-homogeneous  deformation. 

Let  be  "the  coordinates  of  a  particle  referred 
to  a  convected  coordinate  system 


and 


where  r,  0,  and  z  are  the  coordinates  of  the  particle  in  C.3.  The 
Cartesian  coordinates  of  the  particle  in  C.3  are 


x-j  cos  ^  s 


x2  =  sin  ?2  » 


and 


(6.1.1) 
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Let  <J)  be  the  angle  of  twist  per  unit  length  of  the  cylinder. 
The  coordinates  of  the  particle  in  C.l  referred  to  the  Cartesian 
coordinate  system  are 


^1  ~  ^  cos  (^“^g)  » 

X2  =  ^  sin(52-<|>53)  , 

and  X3  =  £3  .  (6.1 .2) 

It  follows  from  equations  (3.2.5)  and  (3.2.6)  that  the  metric 
tensors  in  C.l  and  C.3,  referred  to  the  convected  coordinate  system 
have  covariant  and  contravariant  components 


1  0  0 

0  (S-,)2  -4.(5,  )2 

0  -4(S-,  )2  1+(4>51)2 


1 


0 


0  4>2+(^-)2 

*1 


0 


0 

1 


1  0 

0  (5,)2 

0  0 


0 


1 


(6.1.3) 


< 
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and 


1  0  0 

0  (1/5-,  )2  o 
0  0  1 


(6.1.4) 


Consider  a  particle  at  o'1  in  C.3  as  shown  in  Figure  6.2. 

Let  oxj  be  a  local  Cartesian  coordinate  system  with  origin  o'  and  with 
axes  in  the  directions  of  the  tangents  to  the  convected  coordinate 


Figure  6.2 

Local  Cartesian  Coordinate  System 


curves  passing  through  o'.  If  FiK  is  the  deformation  gradient  tensor 
at  o'  referred  to  the  coordinate  system  ox^ ,  it  may  be  deduced  from 
equation  (3.2.5a)  that  F^F^  are  the  physical  components  of  the  metric 
tensor  G1^  at  point  o'.  These  physical  components  are  given  by 


. 
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G 


=  G 


ij 


^77 


» 


where  the  bars  under  the  repeated  indices  i  and  j  indicate  that  there 
is  no  summation  on  these  indices.  Thus 


1  0  0 

0  1+tS^)2  5^ 
0  5^  1  . 


This  indicates  that  for  the  deformation  defined  by  equations  (6.1.1) 
and  (6.1.2)  each  material  element  in  the  cylinder  undergoes  a  simple 
shear  deformation  with  the  parameter  K,  as  defined  in  Chapter  V,  equal 
to  £-j<J>. 

6.2  Stresses  Prior  to  Yielding  and  in  the  Elastic  Region  After  Yielding 

The  solution  for  the  stresses  required  to  maintain  a  state 
of  torsion  in  a  right  circular  cylinder  of  a  neo-Hookean  material  was 
first  obtained  by  Rivlin  [63],  This  solution  shows  that  a  deformation 
given  by  equations  (6.1.1)  and  (6.1.2)  can  be  maintained  by  surface 
tractions  on  the  ends  of  the  cylinder  without  body  forces  or  other 
surface  tractions . 

Consider  a  right  circular  cylinder  composed  of  a  neo-Hookean 
el astic-perfectly  plastic  material,  subjected  to  the  deformation  de¬ 
scribed  in  section  6.1.  Assume  that  yielding  has  occurred  in  the 
region  a  <_  £-|  <_  a  and  that  for  0  <_  ^  <_  a  the  material  remains  elastic. 


. 

' 

. 

■ 
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that 


Since  K  =  £,<j>  it  follows  from  equations  (5.2.2)  and  (5.2.3) 


for  materials  which  obey  the  von  Mises  and  the  maximum  shear  strain 
energy  yield  conditions  respectively. 

In  the  elastic  region  the  convected  stress  tensor  is  given 

by 


T 


•  • 
1J 


o  •  «  • 

2C  G1J  +  pg1J 


since  the  material  is  assumed  to  be  neo-Hookean  prior  to  yielding. 
From  equations  (6.1.3)  and  (6.1.4)  this  leads  to 


T 


11 


2C+p  , 


t22  =  2C[<t>2  +  1  ?]  +  — 2- 

(S-,)2  (5,) 

23 

t  J  =  2C(j), 


T 


12 


13 


T 


r  (6.2.1) 


and 


=  0  . 


j 


. 
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Substitution  of  the  stresses  (6.2.1)  into  the  equilibrium 
equations  in  cylindrical  polar  coordinates  (Appendix  B)  gives 


H7  -  2C^i  ■ 


0  , 


0  , 


and 


iJL-  = 


0  , 


from  which 


dp(Si)  o 

^57 ' 2C*  5i =  0 


Integration  with  respect  to  ^  yields 


P(CP  =  C<f.2(51)2  +  c1 


vtere  the  integration  constant  c-j  is  determined  from  the  condition 

that  the  radial  stress  is  a  continuous  function  of  the  radius. 

Let  a*  be  the  as  yet  unknown  radial  stress  at  the  elastic-plastic 

boundary  ^  =  a,  and  let  be  the  stress  tensor  at  a  point  referred 

to  the  local  Cartesian  coordinate  system  oxj  at  that  point.  That  is 

a.,  are  the  physical  components  of  the  convected  stress  tensor.  Thus 
1 J 


- 


If  tir  2  ,3  J  2 HO J 

. 
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a- |i  =  2C+p 
=  2C  +  Ct^)2  +  c] 

so  that  a-j  -|  (a)  =  a* 

implies  that  c-j  =  a*  -  2C  -  Ca2c^  . 

In  the  elastic  region  then,  the  stress  solution  in  non-dimensional ized 
form  is 


and 


°11  C  r,r>2  ,a,21i(2  .a*  a33 

tt  k  [(?>  -  (IJ  >  +  T  *  ~T 


a22  .  3C  ,r\2  2  C  ,a,Z 

—  -  -  (¥)  ¥  -  r  W  + 


k  'a  k  * 


23  _  2C  ,r\  , 
(t) 


k  V 


a12  al 3 


=  0  , 


-  (6.2.2) 


where  ip  =  a4> .  The  solution  prior  to  yield  is  obtained  from  (6.2.2) 
by  setting  a  =  a  and  a*  =  0  . 

Let  N  and  M  be  the  total  end  force  and  twisting  moment  re- 
0  0 

spectively  acting  of  the  ends  of  the  cylinder  in  the  elastic  region 


' 


■ 


0  <  Si  £  a.  Thus 


and 


a 

► 

2tt^1°33  d^l 

« 

o 

fa  2 

.  2^/o23d5l 

0 


so  that  substituting  for  a ^  and  °23  from  e<luation  (6.2.2),  performi 
the  integration  and  putting  in  non-dimensional ized  form  gives 


N 


ka 


r  C  ,2/Ou4  ,  0*  fOLst-i 

w[-  2k  *  (I>  +  T  ] 


■exp¬ 


and 


M 

e  _  ttCiLi  /Ou4 
ka3  k  la)  ' 


6.3  Stresses  in  the  Elastic-Plastic  Region 

As  shown  in  section  6.1,  the  deformation  of  each  material 
element  in  the  cylinder  is  equivalent  to  a  simple  shear  deformation. 
Thus  in  the  elastic-plastic  region  the  stress  deviator  tensor 
referred  to  the  local  Cartesian  coordinate  system  is  known  in  tabu¬ 
lated  form  from  the  simple  shear  solution  found  in  Chapter  V. 

The  total  stress  tensor  is  given  by 


a . .  =  a !  .  +  p6 .  •  , 

ij  iJ  iJ 


where  P  =  j  is  the  unknown  mean  normal  stress  which  must  be 


. 

■  :  n 

r  :  n  i  .  '  -  jeJ'rf 
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determined  from  the  equilibrium  equations.  This  yields 


all  (°ii+P)  +  5^  (aira22*  '  0 


(6.3.1) 


and 


=  0 

9^2  9^3 


Integration  of  equation  (6.3.1)  gives 


°ii  +  p  =  - 


-j  [a-j^p)  -  d22(p)]dp  +  c2 


and  since  the  surface  £-j  =  a  is  assumed  to  be  stress  free 


Ca  1 1  ( P )  "  a22^p^dp 


and 


p(C-j)  =  -  crj^^)  + 


a 

1 


p  [a-j  -j  (p)  -  a22 ( P )  * 


1 


Thus  in  the  elastic-plastic  region  the  stress  solution  is  obtained 
from 


a 

1 


’ll 


P  Ha-j  (p)  -  a^2( p)]dp  =  , 


’33 


1 


a22  =  a22~al 1  + 


a 

r  i 


[a-j  i  (p)  -  a22(p)]dp 


a23  *  a23  ‘ 


-  (6.3.2) 


j 


«* 


' 


- 
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The  constant  a*  defined  in  section  6.2  is  obtained  from  the 
above  also  and  is 


a 


”•  la\-\  (p)  -  cr£2(p)]  dp  . 


The  total  end  force  and  twisting  moment  respectively,  acting 
on  the  §nds  of  the  cylinder  are  given  in  non-dimens ionali zed  form  by 


and 


2*(t> 


°33 

k 


^1  /a 
a/ a 


?!  2  ° 


^(-r) 


23 


d(-i) 


(6.3.3) 


(6.3.4) 


6.4  Numerical  Solution 

The  solution  to  this  problem  is  necessarily  a  numerical  one 
since  it  is  based  on  the  previously  obtained  numerical  solution  to  the 
elastic-plastic  simple  shear  problem.  The  latter  solution  gives  the 
stress  deviator  tensor  as  a  function  of  the  parameter  K  and  since  for 
the  torsion  problem  K  =  (—  )ip  the  integrals  in  equations  (6.3.2)  may 

a 

be  evaluated  numerically  using  the  tabulated  values  for  a.j  .  . 

A  computer  program  has  been  written  to  do  the  numerical 
integrations  in  equations  (6.3.2),  (6.3.3),  and  (6.3.4).  The  integrations 
are  preformed  using  the  sixth  order  Gaussian  quadrature  technique  [67] 
described  in  Appendix  D.  The  stress  a.,  is  determined  at  a  finite 

'  J 


' 

to  nottonut  6 

1 1  I 
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si 

number  of  discrete  points  along  a  radius  at  intervals  of  A(— )  equal 
to  0.05.  The  computation  proceeds  from  ^/a  =  1  to  C-j/a  =  a/a  so 
that  at  any  of  the  computation  points  in  the  elastic-plastic  region 
the  integral  in  equation  (6.3.2)  is  given  as  a  sum  of  the  previously 
calculated  integral  plus  the  integral  over  the  last  A(£-|/a)  interval. 

Since  the  solution  for  the  stress  deviator  tensor  which 
appears  in  the  integrand  is  known  only  for  discrete  values  of 
K  =  (C-| /a)ip »  it  is  necessary  to  perform  an  interpolation  to  deter¬ 
mine  the  values  at  the  intermediate  values  of  K  as  required  by  the 
Gaussian  quadrature  routine.  This  interpolation  is  done  using  the 
fifth  order  Newton  forward  formula  discussed  in  Appendix  E. 

The  use  of  the  Gaussian  quadrature  method  for  the  integrations 
in  equations  (6.3.3)  and  (6.3.4)  requires  also  that  interpolations  be 
made  between  the  tabulated  values  of  the  stress  solution  a...  It  is 

*  \J 

not  possible  to  use  Newton's  forward  formula  to  perform  this  interpo¬ 
lation  near  (£-|/a)  =  1  since  the  higher  differences  require  stress 
points  outside  the  actual  cylinder.  This  difficulty  may  be  overcome 
by  using  a  backward  difference  formula  but  the  additional  problem 
arises  that  the  interval  [a/a,l]  may  not  at  times  just  after  yielding 
contain  enough  stress  points  to  allow  the  use  of  a  fifth  order  inter¬ 
polation  formula.  Since  the  stress  solution  is  quite  smooth,  first, 
second,  and  third  order  interpolation  formulae  are  used,  with  the  first 
and  second  order  formulae  being  used  only  when  the  number  of  stress 
points  in  the  interval  [a/a,l]  prohibits  the  use  of  the  third  or 


. 

■ 

■ 
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higher  order  formulae. 

6.5  Discussion 

Some  of  the  numerical  results  obtained  for  both  the  von  Mises 
and  the  maximum  shear  strain  energy  yield  conditions  are  shown  graphically 
in  Figures  6.3  to  6.15.  The  solutions  have  been  found  for  values  of 
k/C  equal  to  0.1,  0.5  and  1.0  and  for  the  equivalent  values  of  W*/C  with 
taking  the  values  0.75  and  1.5. 

The  shapes  of  the  stress/k  versus  r/a  curves  are  of  course 
quite  similar  to  the  shapes  of  the  stress  deviator/k  versus  K  curves 
of  the  simple  shear  solution. 

For  values  of  k/C  <<  1  the  stress  solution  is  that  of  the 
classical  small  strain  el asti c-perfectly  plastic  theory. 

As  the  deformation  parameter  ip  reaches  values  of  two  or 
three  times  that  of  k/C,  the  normal  end  force  and  the  twisting  moment 
approach  constant  values  and  ultimately  they  become  independent  of 
subsequent  twisting  as  does  the  twisting  moment  in  the  classical 
theory  when  almost  the  entire  cylinder  has  yielded. 

It  is  also  observed  that  as  k/C  goes  from  0.1  to  1.0  the 

3 

resulting  change  in  the  value  of  M/ka  for  large  ip  is  only  of  the 

order  of  magnitude  of  0.1  whereas  the  resulting  change  in  the  value 
o 

of  N/ka  for  large  ij;  is  of  the  order  of  magnitude  of  10.  That  is, 
as  in  the  simple  shear  problem  it  is  the  normal  stress  effects  which 
change  radically  as  the  value  of  k/C  increases. 
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Figure  6.3  von  Mises  Yield  Condi tion ,a. ./k  versus 


S-l  03 


STRESS/k 


105 


L 

Q 


Figure  6.4  von  Mises  Yield  Condition,  a-jj/k  versus  j 
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e  6.5  von  Mises  Yield  Condition,  a-./k  versus 
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6.8  von  Mises  Yield  Condition,  a.  ./k  versus  - 
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Figure  6.9  von  Mises  Yield  Condition,  N/ka  versus 
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Figure  6.10  von  Mises  Yield  Condition,  M/ka  versus 
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Figure  6.11  Maximum  Shear  Strain  Energy  Yield  Condition, 
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Figure  6.12  Maximum  Shear  Strain  Energy  Yield  Condition, 
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Figure  6.13  Maximum  Shear  Strain  Energy  Yield  Condition, 

a  I  ./k  versus  ~ 
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Figure  6.14  Maximum  Shear  Strain  Energy  Yield  Condition, 

a../k  versus  -7 
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Figure  6.15  Maximum  Shear  Strain  Energy 
Yield  Condition  a.  ./k  versus  r/a 
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Figure  6.16  Maximum  Shear  Strain  Energy 
Yield  Condition,  o../k  versus  r/a 
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Figure  6.17  Maximum  Shear  Strain  Energy  Yield  Condition, 
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Figure  6.18  Maximum  Shear  Strain  Energy  Yield  Condition, 
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CHAPTER  VII 

EXPANSION  OF  A  THICK  SPHERICAL  SHELL 


7.1  Introduction 

The  problem  of  the  expansion  of  a  spherical  shell  has  received 
considerable  attention  due  to  its  technological  importance.  The  classical 
Hookean  elastic  solution  was  obtained  in  1852  by  Lam£  [68]  and  the  solu¬ 
tion  for  an  hyperelastic  material  was  obtained  in  1950  by  Green  and 
Shield  [69],  The  first  correct  solution  for  the  el asti c-perfectly 
plastic  problem  was  found  by  Reuss  [5].  In  that  analysis  the  total 
deformations  were  assumed  to  be  small.  Hill  [70],  [71]  obtained  a 
solution  to  the  elastic-perfectly  plastic  problem  for  large  total 
deformations  but  retained  the  assumption  that  the  elastic  strains 
are  smal 1 . 

A  solution  is  obtained  in  this  chapter  for  the  stress 
distribution  in  an  expanded  thick  spherical  shell  composed  of  a 
material  whose  elastic  behavior  is  neo-Hookean  and  which  yields 
according  to  any  isotropi c  yield  condition.  Due  to  spherical  sym¬ 
metry  the  stress  at  any  point  in  the  shell  is  a  uniaxial  compression 
-a  superimposed  on  a  hydrostatic  tension  so  that  all  isotropic  yield 
conditions  become  o  =  Y  where  Y  is  the  uniaxial  yield  stress. 

Since  the  deformation  of  the  sphere  is  non-homogeneous , 
residual  stresses  result  from  unloading  from  an  elastic-plastic  state. 
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A  solution  for  the  residual  stresses  is  also  found. 

The  assumption  of  incompressibility  greatly  simplifies  the 
problem.  In  the  elastic-plastic  region  of  the  shell,  the  yield 
condition  and  the  equilibrium  equation  are  sufficient  to  determine 
the  stress  solution  and  it  is  not  necessary  to  consider  the  plastic 
flow  rule  developed  in  Chapter  IV. 

7.2  Prel imi nary 

Referring  to  Figure  7.1  let  (r,e,cj))  be  a  spherical  polar 
coordinate  system  and  ox*  a  Cartesian  coordinate  system. 


x 


Figure  7.1 

Thick  Spherical  Shell  and 
Spherical  Polar  Coordinate  System 


•  _ 
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Four  configurations  of  the  shell,  C.l,  C.2,  C.3,  and  C.4, 
as  defined  in  section  2.1,  are  considered. 

For  a  given  state  of  deformation  C.3,  the  solution  if  no 
yielding  has  occurred,  is  that  found  by  Green  and  Shield  [69].  Since 
the  stresses  are  highest  at  the  inner  surface,  it  is  assumed  that 
yielding  begins  there.  After  the  onset  of  yielding  let  a  be  the  radius 
of  the  elastic-plastic  boundary  which  is  spherical  due  to  symmetry 
and  let  the  radii  of  this  material  surface  in  C.l  and  C.4  be  A  and  a 
respectively.  Similarly  let  R  and  k  be  the  radii  of  the  spherical 
surfaces  in  C.l  and  C.4  which  contain  the  same  material  particles  as 
the  spherical  surface  in  C.3  with  radius  r.  The  inner  and  outer  radii 
of  the  shell  are  then  denoted  by  (R-|  ,R2) »  (r-j,^)  and  (/i-j  ,/l^)  in  con¬ 
figurations  C.l,  C.3,  and  C.4  respectively. 

Define  the  function  Q  by 


From  incompressibility  it  follows  that 


and 


Q(r)  =  [1  + 


R- 


3 

1/3 

] 


(7.2.1) 


The  value  of  Q  at  r  =  ^  is  denoted  by  Q-j  and  is  a  measure  of  the  total 
deformation  of  the  shell  in  C.3. 
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Let  ( 9^3 )  be  coorc|inates  of  a  particle  referred  to 
a  convected  coordinate  system  [48]  defined  by 

*1  =  r  • 

C2  =  e  , 

“  ■*  > 

where  (r,0,<J>)  are  the  coordinates  of  the  particle  in  C.3.  The  co¬ 
variant  and  contravariant  components  of  the  metric  tensor  referred  to 

.  o  »  0  0  0  0 

the  convected  coordinate  system  are  (G..,G1J),  (G.  .,G1J)  ,(g.  and 

vJ  vJ  vj 

•  • 

1 J )  in  configurations  C.l,  C.2,  C.3,  and  C.4  respectively. 

'  J 

Referred  to  the  Cartesian  coordinate  system  the  coordinates 
of  a  particle  in  C.3  are 


x-j  =  C1  sin  cos  C3  » 

sin  sin  £3  > 

X3  ~  ^-j  cos  C2  » 

and  in  C.l  the  coordinates  are 

X]  =  S-|Q  sin  cos  C3  , 


■ 


•• 
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X2  =  sin  C2  sin  , 
=  £,Q  cos  • 


From  equations  (3.2.5)  and  (3.2.6)  and  the  equation 


<10.  =  1  (J _ Q) 

dr  r  q2  » 


(7.2.2) 


the  contravariant  components  of  the  metric  tensor  in  C.l  and  C.3  are 
found  to  be 


[G1j]  = 


0 


0 


1 


2n2 
r  Q 


0 


0 

1 

2  2  2 
r  Q  sin  1 


(7.2.3) 


and 


1  0  0 


[g1J]  ■ 


0 

0 


0 


0 

1 

2  .  2 
r  sin  6 


(7. 2, .4) 


7.3  Elastic  Region 

•  • 

Let  I1"5  be  the  components  of  the  stress  tensor  in  C.3  re¬ 
ferred  to  the  convected  coordinate  system.  The  stress  in  the  elastic 
region  a  <_  r  £  r-j  is  determined  from  the  constitutive  equation 
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•  •  •  •  •  • 
x1J  =  2CG1J  +  pg1J  . 


From  equations  (7.2.3)  and  (7.2.4)  this  gives 


11  _  opn4  , 

T  =  2CQ  +  p  , 

(7 

22  2C  ,  p  .  2q  33 

T  =  o  o  +  9  =  Sin  0T  , 

r  Cr  r 

(7 

and 

12  13  23  n 

(7 

The  equilibrium  equations  (Appendix  C)  are 

22  •  2Q  33  2  11  n 

x  -  r  t  -  r  sin  0x  +  —  t  =0, 

,1  r 

(7 

22  33  22 

x  ?  -  sin0cos0T  +  cot0x  =  0  , 

(7 

and 

33  n 

T  ,3  ‘  0  ' 

(7 

Substitution  for  the  stresses  in  equation  (7.3.4) 

gi  ves 

dp  .  4C  rQ6-2Q3+l i 
dr  r  L  q2  j 

(7 

3.1) 

3.2) 

3.3a,b,c) 

3.4) 

3.5) 

3.6) 

3.7) 


Equations  (7.3.5)  and  (7.3.6)  are  satisfied  identically  due  to  symmetry. 

Using  equation  (7.2.2),  the  first  equilibrium  equation  (7.3.7) 


' 


becomes 


3q  =  4C0-Q3)  . 


integration  of  which  gives 


P  = 


CQ4  +  4CQ  +  c3  , 


where  c^  is  a  constant  of  integration. 

It  is  further  assumed  that  the  outer  surface  of  the  shell 
is  stress  free  so  that  =  0  at  r  = 

Letti ng 


Q2  =  Q(r2) 

it  follows  from  (7.3.1)  that 

2CQ24  +  p(r2)  =  0 


and  from  (7.3.8) 


Thus 


. 


- 


' 
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p  =  -  C(Q4+Q24)  +  c(q-q2)  , 

and  the  mixed  components  of  the  convected  stress  tensor,  which  in  this 
problem  are  identical  to  the  physical  components,  are  after  being 
non-dimensional ized 


and 


TVk'  =  P  (Q4-Q24)  +  F  ’ 


t22/k 


I  _ 


2C  _1_ 
k 1  2 

K  qZ 


P  (q4+q24)  +  p-  (Q-Q2) 


T  3/k'  • 


(7.3.9) 


(7.3.10) 


At  the  elastic-plastic  boundary 


where  Y  =  2k'  ,  k'  being  the  Tresca  yield  stress  in  pure  shear.  Thus 
from  equations  (7.3.9)  and  (7.3.10)  the  equation 

Q36  +  TQ32  -  1  =  0  (7.3.11) 

is  obtained,  where  Q3  =  Q ( a ) .  The  value  of  which  must  lie  in  the 
interval  (0,1)  may  be  found  numerically  from  equation  (7.3.11)  for  a 
given  value  of  k'/C.  Once  is  known,  the  value  of  a  is  found  from 
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equation  (7.2.1)  which  at  r  =  a  may  be  rearranged  to  give 


a 


The  stresses  in  the  elastic  region  and  the  radius  of  the 
elastic-plastic  boundary  can  thus  be  found  for  any  given  state  of 
deformation  as  specified  by  . 


7.4  Elastic-Plastic  Region 

The  yielded  region  of  the  shell  in  C.3  is  given  by 


^  1  r  <  a  , 

and  in  this  region  the  elastic  strains  are  found  using  C.2  as  the 
natural  state.  It  has  been  noted  previously  that  in  general  C.2 

4  « 

corresponds  to  a  non-Eucl idean  space  so  that  the  components  G1J 
cannot  be  found  by  a  transformation  of  the  metric  tensor  from 
Cartesian  coordinates  to  the  convected  coordinates  as  was  done  for 

o  •  o  • 

G1J  and  g1J.  Instead  the  neo-Hookean  constitutive  equation,  the 
equilibrium  equations,  the  incompressibility  of  the  material,  and 
the  plastic  yield  condition  determine  the  components  of  G1J. 

In  the  elastic-plastic  region  of  the  shell  ,  r-j  <_  r  £  a, 
the  stress  is  determined  by  the  yield  condition  and  the  equilibrium 
equation. 

Spherical  symmetry  gives  that 


. 
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T 


12 


=  0  . 


Thus  the  first  equilibrium  equation  (Appendix  C)  becomes 


or 


di 


11 


+  -  t11  -  rr22  -  r  sin20T33  =  0 
dr  r 


dr 


)  =  0  . 


(7.4.1) 


The  stresses  must  satisfy  the  yield  condition 


so  that  equation  (7.4.1)  becomes 


=  4T 
dr  r 


and  after  integration  is 


with 


=  4k' In  r  +  c^ 


+  2k'  . 


The  constant  of  integration  c^  is  determined  from  the  condition 


' 

• 
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that  i  i  be  continuous  at  r  =  a.  The  resulting  stress  solution  is 


t11/k'  =  4  In  (£)  +  (Q34  -  Q24)  +  p-  (Q3  -  Q2)  . 


x22/k'  =  x33/k'  =  T11/k'  +  2  . 


•  • 

The  components  of  the  metric  tensor  G1J  in  C.2  may  now  be 
found.  In  the  yielded  region  the  natural  state  is  C.2  so  that 


a  a  o  o  o  • 

t1j  =  2CG1J  +  pg1J 


(7.4.2) 


12  13  23 

and  therefore  G  ,  G  ,  and  G  are  zero. 

The  condition  that  the  deformation  from  C.2  to  C.3  occur 
without  volume  change  gives 


g/G  =  1 


where  g  =  det(g.  •)  and  G  =  det(G. .)  so  that 

'  J  '  vJ 


01  W3  r4sin20  =  1 


(7.4.3) 


Substitution  of  equation  (7.4.2)  into  the  second  equilibrium 
equation  (7.3.5)  gives 


-J-|£-+  2CG22C0te  -  2CG33sin0cos6  =  0 


- 


. 
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and  since  3p/30  is  zero  from  symmetry 


/sine. 


(7.4.4) 


Thus  from  equations  (7.4.3)  and  (7.4.4)  it  follows  that 


(7.4.5) 


The  yield  condition  and  equations  (7.4.2),  (7.4.4),  and  (7.4.5)  give 


or 


2k' 


(Gll)3/2  +  ^(Gl  1)1/2  .,.0 


(7.4.6) 


This  implies  that  is  constant  throughout  the  yielded  region  and 
comparison  of  equations  (7.3.11)  and  (7.4.6)  shows  that 


G11  =  Q34  .  (7.4.7) 

The  constant  value  of  G1^  determined  from  equation  (7.4.7)  is  hence- 
forth  denoted  by  y. 

A  complete  solution  for  the  stresses  in  the  shell  has  thus 
been  found  for  any  given  state  of  deformation  in  C.3.  It  remains  now 
to  determine  the  residual  stresses  which  result  from  removal  of  the 
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internal  pressure. 

2.5  Residual  Stresses 

It  is  convenient  to  introduce  another  convected  coordinate 
system  defined  in  a  manner  similar  to  the  definition  of  the  coordinates 
Define  the  coordinates  £1  by 

q  =  A, 

V2  = 

and  .  ^3  =  ^ 

where  (a, 0,4)  are  the  coordinates  of  a  particle  in  C.4  referred  to  the 

o  •  «  • 

spherical  polar  coordinate  system.  Furthermore  let  G,1J,  G|1J,  and 
•  • 

be  the  contravariant  components  of  the  metric  tensor  referred  to 
the  convected  coordinates  £!  in  C.l,  C.2,  and  C.4  respectively,  and 
denote  by  t1J  the  contravariant  components  of  the  stress  tensor 
in  C.4  referred  to  this  convected  coordinate  system. 

Since  the  components  of  the  metric  tensor  in  C.2  referred 

i  ^  • 

to  the  £.  coordinate  system  are  known,  the  components  G  J  may  be 


r.ij  _  3E'1  3£'j 
3^  3^ 


found  from 


* 


' 

•:  '\tonsb 
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and  since 


iK-i- 

9£J 


d k 
dr 

0  1 


0  0 
0 


0  0  1 


if  follows  that 


[G,ij]  = 


(— )2y 
^ dr ^  Y 


0 

0 


0 


1 


r2/^ 

0 


0 

0 

1 


r^/y'sin20  . 


The  components  of  g|1J  are 


[5,1J]  - 


1  0 


0 


0  0 


0 


0 


1 


2  2 
A.  sin^e 


The  elastic  deformation  between  C.2  and  C.4  is  isochoric  so 


that 


fl'/G'  =  1  , 


N. 
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where  g'  =  det(gij)  and  G'  =  det(Gjj) 

and  therefore 


(7.5.1) 


Letting  T  =  Vr,  the  condition  of  incompressibility  gives 


TW  ~ 


1 


3  3 

U  -  (*4*-) 


1/3 


(7.5.2) 


Thus  the  components  of  G,1J  are,  after  use  of  equations  (7.4.6)  and 
(7.5.1) 


[G,1J]  = 


JL 

T4 

0 


0 


0 


T  (  4.  I<^ 

*2  (y  c  ) 


0 


0 


0 


I2(y  +  £•) 

o  2 
/i  sin  0 


In  the  region  of  the  shell  in  C.4  which  consists  of  the  ma¬ 
terial  which  yielded  in  C.3,  the  elastic  strains  are  referred  to  the 
natural  state  C.2  so  that 


«  *  •  •  «  • 
t1J  =  2CG,1J  +  pg'1J 


* 
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and  therefore 


and 


t11  =  2C  4  +  p 


.22  _  T  /  ,  k ' \  ,  •  p  _  •  2fi.33 
t  -  2C  ~2  w  +  ~Q~)  +  2  ~  sin  0t 

nr  n. 


(7.5.3) 

(7.5.4) 


The  equilibrium  equations  (7.3.5)  and  (7.3.6),  referred  to 
the  coordinates  are  satisfied  identically  due  to  symmetry,  and  the 
first  equilibrium  equation  gives 


(7.5.5) 


From  the  definition  of  T  and  equation  (7.5.1) 


d  ,1,  4  dT  _  4  (T3- 1 ) 

d n  j4  '  t5  in  -4 


so  that  from  equation  (7.5.5) 


dp.  =  4Cx  (T3-1)2  4k'T2 

tin.  n.  j4  n. 


^  =  4Cy  (H-) 


dT 


T' 


and  thus 
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Integration  gives 


3 

4T  -1 


2k 


■1+T+r 


p  =  Cy(  -j—)  +  1  n[ — — « 

T4  3  (T-l j 


] 


4k1  /2T+i  \  ,  r 

tan  (  )  Cr  • 

sr  /r  b 


(7.5.6) 


The  radial  stress  at  the  inner  surface  of  the  sphere  in  C.4 
is  zero  so  that  from  equations  (7.5.3)  and  (7.5.6)  the  integration 
constant  may  be  determined.  This  results  in  the  following  stress 
solution  for  the  region  ^  <_  n.  <_  a, 


ti  /k.  =  2l(J _ L.)  +  Sx  (1  _  I_) 

r  l/k  k'  M  T  4j  +  k'  lT  T-j ' 


T 


1 


2  0+T+T2)(Trl)2 

+  4  ln  t - 2 - 2*1 

3  O+T^y^d-l)2 


+  —  [tan-1  (^±1)  -  tan’1  (^1)] 


V? 


-/r 


(7.5.7) 


and 


t2  /k'  =  t3,/k'  =  t1  ,/k'  -  (y  +  kl)  (  (7.5.8) 


k'T 


where  T-j  =  T(^-| )  . 

The  value  of  a  is  known  from  the  solution  in  C.3  and  a,  which 
appears  in  T,  is  to  be  determined  from  the  continuity  of  t^  at  n.  =  a. 


- 
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In  the  region  a.  <_> l  <_  ^  the  elastic  strains  in  C.4  are  re¬ 
ferred  to  the  natural  state  C.l  so  that  the  neo-Hookean  constitutive 
equation  in  this  region  is 


tij  =  2CG,lj  +  pg,ij'  . 

•  •  0  0 
The  components  of  G|1J  ,  in  analogy  to  the  components  of  G1J 

are  found  to  be 


[GtlJ]  = 


0 


0 


1 


A.  S 


0  0 


0 

0 

1 


2  2  2 
/i  S  sin^Q 


where  S  =  RM  .  From  the  incompressibility  of  the  material  it  follows 
that 


S  = 


1  + 


A3-*3 


1/3 


The  stresses  are  found  in  a  manner  similar  to  that  used  to 
obtain  the  stresses  in  the  elastic  region  of  C.3.  This  gives 


t \/k‘  = 


(S4-S. 


)  +  £  (S-S,) 


(7.5.9) 
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and 


t  0/k 1  =  t%/k'  = 


2C 


k  1 S 


I 


C_ 

k1 


(s4+s. 


’)  + 


4C 

k' 


(s-s2) 


(7.5.10) 


where  =  S(^)  • 

If  the  value  of  a.  is  known,  equations  (7.5.7),  (7.5.8), 

(7.5.9),  and  (7.5.10)  specify  completely  the  stress  solution.  Con- 
1 1 

tinuity  of  t  at  k  -  a.  gives 


y(~V  -  ~r)  +  4y(^-  -  J-) 

V  V  3  1 

,  2k'  .  (1+T3+T32HTrD2  4k.  .1  2T3+1 

+  In  [ - ^ - yj  +  -  [tan  ( - ) - 


3C 


(l+T1+T12)(T3-l)2 


/3C 


-/T 


1  2T1+1 
tan  ]( — — )] 

-/"T 


-  (S34-S24)  -  4(S3-S2)  =  0  , 


(7.5.11) 


where  =  T(a)  and  Sg  =  S(a).  The  only  unknown  in  equation  (7.5.11) 
is  a  and  it  may  be  found  using  a  numerical  procedure.  The  computer 
program  which  was  written  to  evaluate  the  numerical  results  for  this 
problem  uses  an  interval  halving  technique  [72]  to  solve  equation  (7.5.11) 
for  a. 

7.6  Possibility  of  Further  Yielding  on  Unloading 

It  is  found  that  the  assumption  that  yielding  does  not  occur 
during  unloading  from  C.3  is  not  always  valid.  The  alterations  re¬ 
quired  in  the  residual  stress  solution  when  yielding  does  occur  are 
discussed  in  this  section. 


' 


no  ‘  &,  ■  bn  x  ^  -  ;''v 
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If  it  is  found  using  the  analysis  of  section  7.5  that  the 
yield  condition  is  violated  by  the  residual  stress  solution  then  the 
solution  in  the  region  x^  <_  x  <_  a.  must  be  considered  in  two  parts. 

In  C.4  the  region  in  which  further  yielding  occurs  upon  unloading 
is  assumed  to  be  given  by  • 


1  *  1  b  • 

The  stress  solution  in  the  region 

b  <_  x  <  a 

is  given  by  equations  (7.5.3),  (7.5.4),  and  (7.5.6)  but  the  integration 
constant  c  now  has  a  different  value  which  is  yet  to  be  determined. 

At  x  =  b  the  stresses  satisfy  the  yield  condition 

-  t22  =  2k'  ,  (7.6.1) 


which  gives 


+  T 


kj_ 

Cy 


Ty.  -1=0 


(7.6.2) 


where  =  T ( b ) .  The  value  of  T^  may  be  found  from  equation  (7.6.2) 
using  a  numerical  procedure. 
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Using  equation  (7.5.2)  the  following  equation  is  obtained 
which  determines  the  value  of  b. 


b  = 


a  T 


4‘ 


1/3 


In  the  yielded  region,  ^  <  si  <_  b,  using  the  yield  condition 
(7.6.1)  and  following  an  analysis  similar  to  that  used  to  obtain  the 
stresses  in  the  elastic-plastic  region  of  C.3  gives 


t\/v  =  -  4  In  (^-)  , 


and 


t22/k'  =  t33/k'  =  -4  ln(|-)  -  2 


The  new  value  of  the  integration  constant  c^  in  equation 
(7.5.6)  is  found  from  continuity  of  t^  at  n.  -  b  so  that  in  the  region 
b  <  n  <  a.  the  stress  solution  is 


tli/k'  =  (pr  •  +  4 " 

0+VT4  )(T"1)  /5  -/3 


T? 

■  tan"1 


+  4  In  (-2-) 
*1 


and 


t22/k •  =  t33/k ' 


t1  -|/k 


2Cr  ,  .  k\ 

+  T7"  (Y  + 


2C  Y 

"  k'  t4 


X 


. 


<  J  •  ■  '''  +  ■  1  <'  ■ 
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The  condition  that  the  radial  stress  be  continuous  at  n  =  a 


gi  ves 


•]  +  —  [tan~^  ( — - — )  -  tan” 


-SS 


+  41n  [f-)  -  (S34-S24)  -  4(S3-S2)  =  0  (7.6.3) 


in  which  a  is  the  unknown.  Once  a  is  found  from  equation  (7.6.3)  the 
complete  residual  stress  solution  is  known. 

7.7  Discussion 

A  computer  program  was  written  to  evaluate  numerical  results 
for  the  stress  solutions  obtained  in  this  chapter.  These  results  are 
shown  in  Figures  7.2  to  7.18. 

The  stress  solutions  in  C.3  for  various  values  of  Q-j  and  the 
associated  residual  stress  solutions  are  given  for  shells  with  R^/Ri 
equal  to  2  and  10,  with  values  of  k'/C  equal  to  0.1,  0.5,  and  1.0. 

Figures  7.14  and  7.15  which  give  P/k'  versus  u/R-j  where 
P  and  u  are  the  internal  pressure  and  the  displacement  of  the  inner 
surface  respectively.  All  the  shells  considered  are  seen  to  become 
unstable  before  they  become  fully  plastic.  That  is,  a  point  is  reached 
during  loading  at  which  the  internal  pressure  required  to  maintain  a 


- 
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given  deflection  decreases  with  increasing  deflection.  For  example  a 
shell  with  R^/F^i  =  10  and  k'/C  =  1.0  does  not  become  fully  plastic 
until  u/R-j  equals  8.186  but  becomes  unstable  at  about  u/R-j  =  1.8. 

Thus  the  large  elastic  strains  associated  with  values  of  k'/C  of 
the  order  of  0.1  to  1.0  result  in  significant  changes  in  the  surface 
area  of  the  shells  which  have  a  considerable  weakening  effect. 

The  expansion  of  a  shell  beyond  the  elastic  limit  results 
in  a  strengthening  of  the  shell  after  unloading  in  that  a  higher 
internal  pressure  is  required  to  cause  yielding  in  the  shell  upon 
renewed  loading  because  of  the  residual  stresses  in  C.4.  For  the 
classical  elastic-plastic  problem,  Hill  [70]  has  noted  that  a  shell 
may  not  be  strengthened  by  more  than  a  factor  of  two.  That  is 
P*/PQ  £  2  where  P*  is  the  internal  pressure  in  the  state  C.3  for 
which  the  material  at  the  inner  surface  will  be  just  on  the  point  of 
yielding  if  the  internal  pressure  is  removed,  and  PQ  is  the  internal 
pressure  at  which  yielding  begins  during  the  first  loading  of  the 
shell.  In  the  small  strain  theory  where  superposition  of  stress  solu¬ 
tions  is  possible,  P*  is  also  the  pressure  at  which  yielding  occurs 
during  a  second  loading  of  the  shell.  Reuss  [11]  has  shown  that  this 
maximum  strengthening  can  be  obtained  only  if  R^/Ri  >.  1.701. 

The  solution  found  in  this  chapter  reduces  to  that  from  the 
small  strain  theory.  Using  k'/C  equal  to  0.001  which  is  of  the  same 
order  of  magnitude  as  that  for  most  metals,  the  solution  found  here 
gives  P*/PQ  equal  to  2.00  and  a  critical  shell  thickness  given  by 


. 

' 
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1.81.  For  shells  with  R2/R^  the  order  °f  1  and  with  k'/C  «  1 
the  values  of  a  and  a  are  almost  identical  and  a  small  numerical 
error  in  the  value  of  a  is  likely  to  produce  significant  errors  in 
the  residual  stress  solution.  It  is  believed  that  this  is  the  reason 
for  the  small  discrepancy  in  the  critical  shell  thickness  when 
k'/C  =  0.001. 

In  Figures7.16  to  7.18,  for  shells  with  k'/C  equal  to 
0.1,  0.5,  and  1.0  respectively,  curves  are  drawn  which  give  u Q/ R-j 
and  u^/R-j  as  functions  of  R^/R-j  where  uQ  is  the  inner  displacement 
in  C.3  at  which  yielding  first  occurs  and  u*  is  the  inner  displace¬ 
ment  in  the  state  C.3  at  which  yielding  first  occurs  upon  unloading. 
The  ratio  P*/PQ  is  also  plotted  as  a  function  of  R2/R1 - 

It  is  seen  that  unlike  in  the  small  strain  theory  P*/PQ  is 
not  constant  but  decreases  with  increasing  k'/C  and  for  a  given  k'/C 
decreases  with  increasing  R2/R1  »  reading  a  constant  value  as  R^/R-j 
becomes  of  the  order  of  10  or  more.  Similarly  the  value  of  R^/R-j 
below  which  no  yielding  can  occur  upon  unloading  is  also  dependent 
on  k'/C  becoming  larger  as  k'/C  increases. 
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Figure  7.2  Stress  Solution  in  C.3 
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Figure  7.3  Residual  Stresses 
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Figure  7.4  Stress  Solution  in  C.3 
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Figure  7.5  Residual  Stresses 
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Figure  7.6  Stress  Solution  in  C.3 
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Figure  7.7  Residual  Stresses 
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Figure  7.9  Residual  Stresses 


20 


152 


,  VSS3H1S 


Figure  7.10  Stress  Solution  in  C.3 
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Figure  7.11  Residual  Stresses 
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Figure  7.12  Stress  Solution  in  C.3 
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Figure  7.13  RESIDUAL  STRESSES 
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Figure  7.14  Internal  Pressure/k'  versus  Inner  Displacement/!^ 
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Figure  7.15  Internal  Pressure/k'  versus  Inner  Displacement/^ 
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Figure  7.16  uq/R^  ,  u*/R.|  and  P*/PQ  versus  F^/R-j 
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Figure  7.17  uq/R-j  ,  u*/R^ ,  and  P*/PQ  versus  R^/R-j 
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Figure  7.18  uq/R-j  s  u*/R-|  9  and  P*/PQ  versus  R^/R-j 
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CHAPTER  VIII 

CONCLUDING  REMARKS 


As  indicated  in  Chapter  I,  the  motivation  for  the  development 
of  a  theory  for  elastic-plastic  flow  with  finite  elastic  strains  is 
twofold.  Firstly  there  is  the  possibility  of  finite  elastic  volume 
change  occurring  in  certain  explosive  metal  forming  processes  and  the 
work  of  Lee  [25],  [26],  [27]  is  concerned  in  part  with  this  phenomenon. 

Some  of  the  work  in  this  thesis  is  also  applicable  to  this  problem  since 
although  the  assumption  of  elastic  incompressibility  is  made  in  Chapters 
IV  to  VII,  the  kinematics  in  Chapter  II  and  the  thermo-elastic  considerations 
for  an  elastic-plastic  material  in  section  3.7  are  not  so  restricted  and 
may  be  used  as  a  starting  point  of  a  general  theory  in  which  elastic 
incompressibility  is  not  assumed.  Furthermore  the  development  of  the 
flow  rule  (4.2.4)  may  be  generalized  [76]  to  include  the  possibility  of 
finite  elastic  volume  change. 

Even  though  the  development  of  such  a  theory  is  possible  it 
is  felt  that  the  solution  of  actual  boundary  value  problems  using  such 
a  theory  may  be  intractable  since  it  is  in  general  extremely  difficult 
to  obtain  solutions  in  the  theory  of  finite  strain  elasticity  for  com¬ 
pressible  materials. 

Two  problems  of  interest  for  which  such  a  theory  may  be  re¬ 
quired  are  the  propagation  of  one  dimensional  elastic-plastic  waves  in 
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a  plate  due  to  an  explosive  charge  on  one  of  the  plate  surfaces  and  the 
explosive  expansion  of  a  spherical  cavity  in  an  infinite  elastic-plastic 
medium.  Lee  [26],  [27]  has  discussed  the  first  problem  and  the  latter 
has  been  considered  by  Hunter  and  Crozier  [77]  as  well  as  in  a  com¬ 
prehensive  review  by  Hopkins  [78].  Only  small  elastic  strains  are 
considered  in  these  two  works. 

The  development  of  a  more  general  theory  for  elastic-plastic 
materials  is  also  the  result  of  indications  that  some  elastomers  and 
polymers,  for  certain  ranges  of  temperature,  exhibit  permanent  set 
after  finite  elastic  shear  strains  with  little  if  any  elastic  volume 
change.  This  is  in  contrast  to  the  previously  mentioned  class  of 
problems  in  which  the  reverse  situation  occurs,  that  is  elastic-plastic 
flow  with  finite  elastic  volume  change  and  small  elastic  shear  strains. 
The  possibility  exists  that  an  elastic-plastic  theory  in  which  finite 
elastic  shear  strains  are  considered  may  provide  a  useful  approximation 
to  the  behavior  of  such  materials  and  the  work  of  Chapters  IV  to  VII, 
in  which  elastic  incompressibility  is  assumed,  has  been  done  with  these 
materials  in  mind. 

The  assumption  is  made  in  Chapter  IV  that  the  materials  con¬ 
sidered  are  rate  independent  and  satisfy  Drucker's  postulate.  At  pre¬ 
sent  there  is  no  evidence  to  indicate  that  for  inelastic  deformation 
the  elastomers  mentioned  necessarily  belong  to  this  class  of  materials 
and  although  in  the  limited  testing  cited  in  Chapter  I  there  was  no 
indication  that  these  materials  were  not  in  this  class  the  question  re- 
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mains  the  subject  of  further  experimental  work.  Furthermore  the  two 
yield  conditions  discussed  in  Chapter  IV  are  not  based  on  experimental 
evidence  and  also  remain  to  be  investigated  experimentally.  Temperature 
is  known  to  have  a  significant  effect  on  the  properties  of  these  ma¬ 
terials  and  any  testing  must  certainly  take  this  into  account  especially 
with  regard  to  the  assumption  of  rate  independence.  It  is  concluded 
therefore  that  considerable  experimental  work  is  required  before  the 
work  of  Chapters  IV  to  VII  may  be  linked  with  the  physical  behavior  of 
any  real  materials. 

Although  the  theory  presented  here  is  limited  to  perfect 
plasticity  there  is  no  difficulty  in  principle  in  including  work¬ 
hardening  in  the  development.  There  is  no  experimental  evidence  how¬ 
ever  to  indicate  that  the  work-hardening  of  the  elastomers  mentioned 
is  or  is  not  significant. 

The  neo-Hookean  material  has  been  chosen  as  the  model  for 
elastic  behavior  in  Chapters  IV  to  VII  because  it  provides  a  reason¬ 
able  approximation  for  the  elastic  behavior  of  some  elastomers  and 
because  it  simplifies  the  mathematical  analysis.  The  inclusion  in 
the  theory  of  more  elaborate  models  for  incompressible  elastic  be¬ 
havior  such  as  the  Mooney  material,  although  likely  to  result  in  in¬ 
creased  difficulty  in  solving  actual  problems,  presents  no  difficulty 
in  principle. 

Elastic*  solutions  are  known  for  a  number  of  deformations 

*Here  the  term  elastic  refers  to  the  general  finite  strain  theory 
and  not  the  classical  Hookean  theory. 
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which  are  included  in  the  five  families  of  deformations  [73]  which  are 
possible  in  any  isotropic  incompressible  elastic  material**  under  the 
action  of  surface  tractions  only.  Among  these  are  the  extension,  in¬ 
flation,  and  torsion  (or  shear)  of  a  hollow  cylinder,  and  the  bending 
of  a  cuboid.  An  elastic-plastic  solution  to  a  problem  cannot  be  found 
unless  the  elastic  solution  is  known  and  these  existing  elastic  solu¬ 
tions  provide  a  number  of  problems  which  may  be  solvable  using  a  finite 
strain  elastic-plastic  theory  such  as  that  developed  here. 


**For  a  specific  elastic  material  there  may  exist  deformations  which 
can  be  maintained  by  surface  tractions  alone  but  which  are  not  in¬ 
cluded  in  these  five  families. 
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APPENDIX  A 

A  THEOREM  USED  IN  SECTION  4.4 

The  following  theorem  which  is  used  to  simplify  the  plastic 
flow  rule  developed  in  Chapter  IV  is  now  proved. 

Theorem:  If  A,  B,  and  C  are  square  symmetric  matrices  with  B  positive 
definite,  and  the  principal  directions  of  A  and  B  coincide,  and  if 

A  =  BCB'1  +  B_1CB  (A. 1 ) 

then  the  principal  directions  of  the  matrix  C  coincide  with  those  of 
and  B. 

Proof: 

Let  x  (  a  column  vector)  be  in  the  direction  of  any  one  of 
the  eigenvectors  of  A,  and  therefore  by  hypothesis  it  is  in  the  di¬ 
rection  of  one  of  the  eigenvectors  of  B„  Thus  there  exist  scalars 
A  and  A  such  that 


Ax  =  Ax 


Bx  =  Ax 


(A. 2) 

(A. 3) 

(A. 4) 


and 
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From  equations  (A.l),  (A. 2),  (A. 3),  and  (A. 4) 


=  (BCB'1  +B‘1CB)x 
=  (1  BC  +  AB_1C)x 

=  ({  S  +  XB_1)Cx  .  (A. 5) 


There  exists  a  scalar  a  and  a  vector  y  such  that 


where 


Cx  =  ax  +  y 

'N/  >V  ^ 

/x  =  0  .  (A. 6) 


Thus  equation  (A. 5),. becomes 


Ax  =  (-j-  B  +  AB  ^)(ax  +  y) 

=  (f  §X  +  Aa  B_1x)+(1  By  +  Xf\) 


2ax  +  (1  By  +  XB_1y) 

~  A  ~~  ~  ~ 


(A. 7) 


' 


173 


Premultiplying  equation  (A. 7)  by  and  using  the  condition 
(A. 6)  gives 

/§*  +  ^  -  o  . 

and  since  B  is  positive  definite  y  must  be  the  null  vector  so  that 

Cx  =  ax  . 

'S/'V  »v 

That  is,  x  is  an  eigenvector  of  C.  Consequently  each  of  the  eigen¬ 
vectors  of  A  or  B  has  the  same  direction  as  one  of  the  eigenvectors 
of  C  so  that  A,  B,  and  C  are  coaxial. 
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APPENDIX  B 

EQUILIBRIUM  EQUATIONS  IN  CYLINDRICAL  POLAR  COORDINATES 


Consider  the  cylindrical  polar  coordinate  system  shown  in 
Figure  B. 1 . 


Figure  B.l 

Cylindrical  Polar  Coordinate  System 


Let  the  quantities  £ .  be  defined  by 


_ 
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r  =  , 

e  =  > 

and  z  =  £3  . 

An  element  of  length  ds  is  given  by 

2  2  2  2  2 
ds  =  dr  +  r  d0  +  dz  , 


so  that  the  metric  tensor  referred  to  the  cylindrical  polar  coordi¬ 
nates  E.  has  components 

11  ,  _  1  _  2  33  _  n 

9  "  9n  "  1  »  922  "  22  r  *  g33  “  9  * 

9 

and  all  other  components  are  zero. 

The  Christoffel  symbols 


ININ  ,  J 


are  thus 


r2  -  r2  -  1  r1 
1 12  "  X21  r  *  1 22 


=  -  r  , 


and  all  other  are  equal  to  zero. 


(B.l) 


n  .  ■  ,  .  f:  r 


*"  N  i 


e '  ££e  *  1  ":sr  ’  ss*  •  f  ’  rr®  “  ® 
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Cauchy's  first  equation  of  motion  (3.3.3)  reduces  to  the 
equilibrium  equation 


t’J  -  0 

5  I 


for  a  quasi-static  process  in  the  absence  of  body  forces.  That  is 


Jj  +  r1  TkJ  +  rJ  T1k  =  0 
1  ,1  Jki  +  1  ki  1  u 


so  that  using  equations  (B.l)  the  equilibrium  equations  in  cylindrical 
polar  coordinates  are 


11  .  21  ,  31  .  1  11  22  n 

t  t+t  0  +  t  +  -  ri  =0, 

,1  j2  j  3  r 


12  '22  .  32  ,  3  12  n 

T  ,1  +  T  ,2  +  T  ,3  +  7T  =  0  > 


13  .  23  ,33  .  1  13  n 

T  ,1  +  T  ,2  +  T  ,3  +  7T  =  0  • 
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APPENDIX  C 

EQUILIBRIUM  EQUATIONS  IN  SPHERICAL  POLAR  COORDINATES 
Consider  the  spherical  polar  coordinate  system  in  Figure  C.l. 


Figure  C.l 

Spherical  Polar  Coordinate  System 


Let  the  quantities  E.  be  defined  by 


r 


> 
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6  = 


» 


and 


An  element  of  length  ds  is  given  by 

ds2  =  dr2  +  r2de2  +  r2sin20d4>2 

so  that  referred  to  the  spherical  polar  coordinates  £ .  the  metric 
tensor  has  components 


11  12  12.2 
9  9-| -|  =  1  »  922  “22  r  *  g33  =  ~"33  =  r  sin  6 

9  9 


with  all  other  components  zero, 

* 

The  Christoffel  symbols  are  thus  found  to  be 


v\z  =  -  r  ,  r 


33  rsin  9  »  r!2  r  ’ 


33 


3  1  3 

sin6cos0,  =  -  and  r ^  =  cot6 


with  all  others  being  zero. 

Thus  the  equilibrium  equation 


. 


becomes 


T 


,  referred  to  the  spherical  polar  coordinates  £.  , 


1  +  T21  ?  +  T31  ~  +  cot0T21  +  7T11  -  rr22  -  rsin20i33 
» *  3  r 


12  a.  22  ,  32  ,  4  12  ,  22  .  -  .  33  n 

1  +  t  o  +  T  ^  +  —  t  +  c°t0T  -  sin0cos0i  =  0 


13  ,  23  ,  33  ,  4  13  ,  0.a  23  n 

T  -j  +  t  2  +  t  3  +  Ft  +  0  • 


and 
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APPENDIX  D 

GAUSSIAN  QUADRATURE  INTEGRATION 

The  Gaussian  quadrature  method  of  numerical  integration, 
which  is  known  for  its  simplicity  and  accuracy,  is  discussed  briefly 
in  this  appendix. 

As  a  preliminary  however,  the  Lagrangian  interpolation  poly¬ 
nomial  [74]  for  unequally  spaced  data  must  be  considered. 

If  data  y^,  y2»  . . . ,  yn  are  known  at  points  x-j ,  x2,  ...»  xn 
which  are  not  necessarily  equally  spaced  then  a  polynomial  <j>  -|(x) 

which  is  of  degree  n-1  may  be  fitted  through  these  points.  Such  a 
polynomial  is  the  Lagrangian  interpolation  polynomial 


(x-x2)(x-x3)  ...  (x-xn) 
Vl^  "  "(x1-x2)(x1-x3)  ...  (x1-xn)  yl 

(x-x1)(x-x3)  ...  (x-xn) 

(x2-x-| )  (x2-x3)  .  • .  ( ^  ^ 


+  ... 

(x-x1)(x-x2)  ...  (x-xn i ) 


n 


Vl(x)  =  ^  Li(x)  *1 


or 


(D.l) 


■ 


' 


181 


where  the  Lagrangian  coefficients  L. (x)  may  be  written  as 


Li  (x) 


6(x) 

(x-x- )3' (x.) 


and 


3 ( x )  =  (x-x1)(x-x2)(x-x3)  ...  (x-xn)  . 


Since  an  integral 


I 


b 

* 

f(s)  ds 
a 


may  always  be  put  in  the  form 


+1 

* 

f(x)  dx 


by  the  change  of  variable 


(D.2) 


(D.3) 


2s  -  (a+b) 
(b-a)  ’ 


only  the  integral  (D.3)  is  considered  in  the  following. 

The  question  arises,  how  may  the  integral  I  be  put  in  the 


form 


+1 

» 

f  (x)dx  =  w^x-j)  +  w  2y(x2)  +  ...  +  wny(xn)  (D.4) 


-1 


. 


' 


•  „  j  if  H  I  |  H  fill  irnoM 
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so  that  it  is  exact  for  f(x)  as  high  a  degree  polynomial  as  possible. 

For  n  points  x-j ,  x^,  ...  ,  xn#  a  polynomial  <|>  -|(x)  of 

degree  n-1  determined  by  equation  (D.l)  may  be  fitted  through  the 
corresponding  ordinate  points.  If  the  integrand  f(x)  is  a  polynomial 
of  degree  n-1  or  less  then  4>  i (x)  and  f(x)  are  identical.  However  if 
f(x)  is  a  polynomial  of  degree  2n-l  or  less  it  may  be  written  as 


f(x)  =  V|(x)  +  Vl^  »  (D-5) 


where  $(x)  is  a  polynomial  of  degree  n  as  given  by  (D.2)  and  zn_-|(x) 
is  a  polynomial  of  degree  n-1.  Integration  of  (D.5)  gives 


+1 

» 

f(x)dx 

-1 


n 

v'V +  Vi  > 


where 


+1 

L.(x)  dx 

« 

-1 


(D.6) 


and 


+1 

» 

B(x)Z  i  (x)  dx  . 
-1  n'' 


The  term  Rn  -|  is  the  error  which  results  from  representing  the 
function  f(x)  by  the  polynomial  <j>  -j(x)  .  A  sufficient  condition 
that  Rn  i  be  zero  is  that  $(x)  be  a  Legendre  polynomial  of  degree  n 
since  the  Legendre  polynomials,  which  may  be  found  from 


'  %  m  z\ 


, 


■ 
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n 


Pn(x>  =  (x2_1)n  • 

n  2nn !  dxn 


have  the  property  that 


+1 

»  4 
P  (x)  xmdx  =  0 

-1 


for  all  m  <  n  . 

If  G(x)  is  chosen  equal  to  P  (x)  then  the  values  of  x-j ,  Xg» 
xn  are  known  since  they  are  the  roots  of  the  polynomial  P  (x). 

If  data  are  chosen  at  these  n  points  and  the  weighting  coefficients 
w.  are  found  from  equation  (D.6),  the  integration  formula  (D.4) 
will  be  exact  if  f(x)  is  a  polynomial  of  degree  2n-l  or  less.  This 
partly  explains  the  accuracy  of  the  method  since  most  other  integration 
formulae  using  n  data  points  are  exact  only  if  the  integrands  are 
polynomials  of  degree  n-1  or  less. 

The  roots  x.  of  the  Legendre  polynomials  and  the  weighting 
factors  w.  are  readily  available  in  tabulated  form  [75]. 


. 
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APPENDIX  E 

NEWTON'S  INTERPOLATION  FORMULAE 

Let  y  ,  y-, ,  . . . ,  y  be  the  values  of  the  function  f(x)  at 
J  o  l  n 

n+1  points  xQ,  x-j . xn  which  are  equally  spaced  by  a  distance  h. 

It  is  possible  to  approximate  the  function  f(x)  by  a  polynomial  4>n(x) 
of  degree  n  so  that  values  of  f(x)  at  intermediate  values  of  x  may  be 
estimated.  This  is  the  basis  of  Newton's  interpolation  formulae. 

The  derivation  of  these  formulae  is  simplified  by  the  use 
of  three  operators,  E  the  shift  operator,  A  the  forward-difference 
operator,  and  V  the  backward-difference  operator.  These  are  defined 
by 


E  f(x)  =  f(x+h)  , 

Af(x)  =  f(x+h)  -  f(x)  , 

and  Vf (x)  =  f(x)  -  f(x-h)  . 

Integer  powers  of  the  operators  A  and  V  are  defined  by 


Amf(x)  =  Am_1 f (x+h)  -  Am_1f(x)  , 


■ 


185 


and 


Vmf(x)  =  Vm_1f(x)  -  vm_ 1 f ( x-h )  , 


where  m  is  an  integer. 

Furthermore  integer  powers  of  E  are  defined  by 

Emf (x)  =  f (x+mh) 

and  extending  this  to  non-integer  powers  gives 

Eaf(x)  =  f (x+ah) 

where  a  is  any  real  positive  number. 


Two  important  relations  between  the  operators  E,  A,  and  V 


are 


A  =  (E-l ) 


(E.l) 


V  =  ( 1-E"1 )  . 


and 


These  are  easily  verified  using  the  definitions  of  the  operators  E, 
A,  and  V. 


Consider  the  following  expansion  based  on  equation  (E.l) 


Ea  =  (1+A)a 


I 
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=  1  +  c»A  +  ^  +  01^ct~^a~2^  A3  +  ...  (E.2) 

If  a  is  equal  to  an  integer  value  n,  the  above  series  is  finite  with 
n+1  terms.  Applying  the  operator  (E.2)  at  the  point  y(xQ)  gives 

y(xQ+ah)  =[1  +  aA  +  a--|— —  A^  +  ...  + 

0‘(0‘-'^''ni(g'-n+1)  A"]  yl'xo)  +  R(x)  (E'3) 


where 


R(x)  =  +  •••]  y(xj  . 


Denoting  by  ^(x)  the  first  term  on  the  right  hand  side  of  equation 
(E.3)  gives 


y(xQ+ah)  *  <|>n(x)  +  R(x)  . 

If  a=n,  the  term  R(x)  is  zero  since  the  series  (E.2)  termi¬ 
nates  after  n  terms  so  that  the  expression 

y(xQ+ah)  =  <{>n(x) 


where 


' 


*n(x) 

is  exact, 
fo rmulae. 


gives 

where 

*n(x) 


=  [1  +  aA  +  A2  +  ...  +  An]y(x0)  (E.4) 

Equation  (E.4)  is  called  Newton's  forward  interpolation 
similar  expansion  of 

E01  =  (l-v)"a 
y(xn+ah)  =  <f>n  (x)  +  R(x) 


[1  +  av  +  v2  +  ...  +  0‘(a+1)--;(c‘+n-l)  Vn]y(x  ) .  (E.5) 

Cm  •  I  I  • 


Equation  (E.5)  is  called  Newton's  backward  interpolation  formula. 


